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Abstract 

We calculate the kaon B-parameter, Bk, to next-to-leading order in staggered chiral perturbation 
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' theory. We find expressions for partially quenched QCD with three sea quarks, quenched QCD, 
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£Nj ■ and full QCD with m u = uid / m s . We extend the usual power counting to include the effects of 

, using perturbative (rather than non-perturbative) matching factors. Taste breaking enters through 
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the 0(a 2 ) terms in the effective action, through 0(a 2 ) terms from the discretization of operators, 
and through the truncation of matching factors. These effects cause mixing with several additional 
operators, complicating the chiral and continuum extrapolations. In addition to the staggered 
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expressions, we present Bk at next-to-leading order in continuum PQxPT for Nj = 3 sea quarks 
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with m u = ma ^ m s . 
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I. INTRODUCTION 



Experimental measurements of CP violation can be used to extract information about the 
CKM matrix. In particular, the size of indirect CP violation in the neutral kaon system, 6k, 
combined with theoretical input, places an important constraint on the apex of the CKM 
unitarity triangle . Because 6k is well-known experimentally [3|, the dominant source of 
error in this procedure is the uncertainty in the lattice determination of the nonperturbative 
constant Bk, which parameterizes the K° — K° matrix element. Because new physics would 
likely produce unitarity violation in the CKM matrix and additional CP-violating phases, a 
precise determination of Bk will help to constrain physics beyond the standard model. 

Promising calculations with partially quenched staggered fermions are in progress 
Because staggered fermions are computationally cheaper than other standard discretiza- 
tions, they allow simulations with the lightest dynamical quark masses currently available. 
Unfortunately, staggered fermions come with their own source of error - taste symmetry 
breaking. Each lattice fermion flavor comes in four tastes, which are degenerate in the con- 
tinuum. The nonzero lattice spacing, a, breaks the continuum taste symmetry at 0(a 2 ), and 
the resulting discretization errors are numerically significant at present lattice spacings Q. 
Thus one must use the chiral perturbation theory functional forms for staggered fermions in 
order to correctly perform the combined continuum and chiral extrapolations incorporating 
taste violations jaSQj- It is clear, then, that the potential improvement in precision of the 
lattice determination of Bk is limited without the availability of the appropriate staggered 
%PT expression. 

In this paper we calculate Bk to next-to-leading order (NLO) in an extended version 
of staggered chiral perturbation theory (S^PT) that incorporates perturbative operator 
mixing. In the NLO S%PT expressions for quantities calculated thus far, e.g fx and fjj, 
taste-breaking primarily enters through additive corrections to the tree- level masses of pions 
inside loops [a U . The calculation of Bk in S%PT provides a more complicated example 
of taste-breaking in which new operators in the staggered chiral effective theory significantly 
change the result from that in continuum xPT, even at NLO. The bulk of our new work 
is in correctly enumerating all possible operators which contribute. It turns out that the 



1 Calculations of Bk with dynamical domain wall fermions are also underway 
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NLO expression is in terms of 37 low-energy constants (21 for a single lattice spacing), 5 of 
which are already known from fitting other staggered lattice data. Because fitting so many 
parameters requires a large amount of data, we show how to extract additional coefficients 
using simple kaon matrix elements. Nevertheless, performing the continuum and chiral 
extrapolations of Bk will be challenging. 

We aim to keep the body of this paper accessible to those wishing to use our NLO S%PT 
expressions for calculations of Bk on the lattice. Our paper is therefore organized as follows. 
In Section |n] we define Bk in QCD, and discuss modifications to the quark-level operator 
necessary for its calculation with staggered quarks. We summarize the main results of our 
Bk operator enumeration in Section 11111 In Section II VI we calculate Bk at 1-loop for a 
partially quenched (PQ) theory with three dynamical quark flavors (each with four tastes). 
We then add the NLO analytic terms and give results for Bk to NLO in quenched, partially 
quenched, and full QCD in Section |Vj Finally, in Section EU we explicitly show, for a single 
lattice spacing, how to determine six of the undetermined NLO coefficients, as well as how to 
use the S%PT fit results to extract the renormalization group invariant quantity Bk in the 
continuum. We conclude in Section IVH1 As a corollary to our staggered result, we present 
the expression for Bk at NLO in continuum PQxPT with 2+1 sea quarks (m u = ^ m s ) 
in Appendix [XI since it is not available in the literature. Appendix [B] briefly reviews the 
essentials of staggered chiral perturbation theory necessary for understanding the details of 
our NLO calculation and operator enumeration. Appendix ICl enumerates in detail all of the 
operators which contribute to Bk at NLO. 

II. B K WITH STAGGERED QUARKS 

In this section we relate Bk in QCD to matrix elements of a continuum staggered theory 
in which one has taken the fourth-root of the quark determinant but there are still four 
tastes per valence quark flavor. 

The kaon 5-parameter is defined as a ratio of matrix elements: 

M K = (K^\O k \K°) = B K M vac , (1) 

where Ok is a weak operator and Ai vac is the result for A4k in the vacuum saturation 
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approximation: 

O k = [^7m( 1 + 75)4][^7m( 1 + 75)4] (2) 
-Mvac = ^(^|[^(1 + 75)4]|0)(0|[^7m(1 + 75)4]|/O- (3) 

Note that there are separate summations over color indices a and b. Thus the matrix element 
M.k receives contributions from two different quark- level contractions, one of which produces 
a single color loop and the other which has two color loops. The expression in Eq. {I} can 
be simplified using the fact that -M V ac is related to the square of the kaon decay constant: 

M K = % -B K m\f K , (4) 

where we use the normalization that fx ~ 156 MeV. 

In order to calculate B>k with staggered fermions we must introduce the taste degree of 
freedom, both in the operators (Ok) and in the states (K° and K°). We choose the external 
staggered kaons to be taste P, by which we mean that the lattice meson operator contains 
the pseudoscalar taste matrix £ 5 . Since this is the lattice Goldstone taste, its correlation 
functions satisfy U(1)a Ward identities, so S%PT expressions for its mass, decay constant, 
and other physical quantities are simpler than those for other tastes of PGBs. In addition, 
it is a local kaon on the lattice, and therefore relatively simple to implement numerically. 
Because the external kaons are taste P, the weak operator should also be taste P: 

OT ve = [<7m(1 + 7s) ® [5F(7a*(1 + 7s) ® Zs)d b ] (5) 

However, when this operator is Fierz-transformed, not only do its color indices change, but 
it mixes with other tastes. To remedy this, we introduce two types of valence quarks, 1 and 
2, into the B>k operator: 

Of" ered = 2|p:(7^(l+75)®e5)rfla]p2;(7M(l + 75)®e5)rf2 fe ] 

+ |si;(7 M (l + 75)®6)dl6]p2^(7^(l + 7B)®6)d2a]J (6) 
and take the matrix element between two types of kaons: 

M stag 9 ered = ^ ^ | pstaggered | R 0^ ^ (?) 

where K® is an sldl meson and K\ is an s2<i2 meson. The extra valence quarks require 
explicit inclusion of both color contractions in Q s ^ a 99 ered ^ w ] 1 i] e the overall factor of 2 ensures 



that A4k and J\4 s ^ a " ered have the same total number of contractions. When Q s £ a 99 ered j g 
Fierz-transformed, it has a new flavor structure so it does not contribute to j\yi s ^ a " ered . Thus 
^staggered canno |- rece j ve contributions from incorrect tastes. We choose that the two sets 
of valence quarks have equal masses, i.e. m^i = = fn x and m sl = m s2 = m y , in order 
to simplify future expressions. 

The continuum staggered theory differs from QCD in that it has four copies of each 
valence quark, so M.k an d j\zl K a " ered are not precisely equal. Nevertheless they are simply 
related by overall normalization factors. Specifically, 

M K = ±-Mf" ered (8) 

(U\s llxl sd\K°) = W^(0|s( TM 75®6)d|^), (9) 



where Nt — 4 is the number of tastes per valence quark flavor [12j|. We emphasize that these 
expressions assume that one has already taken the fourth-root of the quark determinant and 



the continuum limit, and that the vDet procedure is valid. We now re-express Bk in terms 
of the staggered matrix elements calculated on the lattice: 

B K = Nt P K P . (10) 



Because the various factors of iV t cancel, Bk in QCD is egua/ to the value of the analogous 
ratio of staggered matrix elements. 

Current calculations of Bk with staggered quarks use the partially quenched approxi- 
mation. The quark content the corresponding PQ theory is quite large. We have already 
seen that staggering requires two sets of d and s valence quarks. Partial quenching adds 
two corresponding sets of ghost quarks as well as three sea quarks, resulting in eleven total 
quark flavors, each of which comes in four tastes. In order to make this completely clear, 
we show the explicit form of the quark mass matrix: 

M = diag{m u I , m^J, m s I, m x I, m y I, m x I, m y I, m x I, m y I, m x I, m y I}, (11) 

sea valence 1 valence 2 ghost 1 ghost 2 

where / is the 4x4 identity matrix. Thus the generalization of Bk to PQ staggered quarks 
is relatively straightforward conceptually, although nontrivial to implement on the lattice. 
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III. GENERALIZED STAGGERED CHIRAL PERTURBATION THEORY FOR 

B K 



The previous section described how to calculate Bx with staggered quarks in the contin- 
uum, but our expression for Bx in S%PT must describe Bx at a 7^ if it is to be used for 
continuum and chiral extrapolations of lattice data. It is therefore necessary to discuss in 
some detail how Bk is actually calculated on the lattice. 

The Bx matrix element, j\y{ s J c a99ered ; receives a contribution from the lattice version of 
Qstaggered^ &g ag f rom Q^gj- lattice operators that are in the same representation of the 
symmetry group that maps a hypercube onto itself 

stag 9 ered,cont = ^aggeredM + « [t ast g p opg 1 + ^ [other taste P S.] 

47T 47T 

+a 2 [all taste ops.] + a 2 [all taste ops.] + . . . , (12) 

where a is the strong coupling constant. The 1-loop perturbative matching coefficients 
between Q s ^ a 99 ered j n the continuum and four-fermion lattice operators are known, and are 

. However, the 2-loop matching coefficients have not 
been determined, so, in order to remain conservative, we consider them to be of order unity 
times a 2 without any factors of 4ir. Current numerical staggered calculations are in fact of 
the following matrix element: 

(K^ p \0^ loop (taste P)\K° 2p ) = Mat, (13) 

where the subscript "1 — loop" and the argument "taste P" indicate that one includes all stag- 
gered lattice operators with taste P that mix with the latticized Bx operator at (9(a/47r), 
i.e. those in the second term on the RHS of Eq. ()12j) . using the appropriate matching co- 
efficients. However, one neglects wrong-taste and higher-order perturbative mixing (terms 
three and four), as well as all operators which arise through discretization effects (term five), 
in Eq. (fT2*|) . Although the expression in Eq. |TH|) differs from the continuum matrix element, 
it reduces to the desired quantity in the continuum limit. Generically, 

Mat = M ont + -rM' + a 2 M" + a 2 M"' + . . . , (14) 

where A^cont is the desired continuum result. The matrix element Ai' comes from neglecting 
taste- violating 1-loop operator mixing, while M" comes from neglecting 2-loop operator 
mixing. Both taste-breaking and taste-conserving discretization errors generate M!" ' . 
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We are now in a position to determine the appropriate power-counting scheme for calcu- 
lating Bk at next-to-leading order in S^PT. Clearly it must incorporate a 2 , a/47r, and a 2 . 
Continuum x?T is a low-energy expansion in both the pseudo-Goldstone boson (PGB) mo- 
mentum and the quark masses; it assumes that Fpg^/A 2 ~ ^pgb/^x' wnere m PGB ^ m i 
and A x pa 4irf n is the %PT scale. However, the numerical values of m q , a 2 and a all depend 
on the particular parameters of a given lattice simulation. Current PQ staggered lattice sim- 



ulations [7] use a range of PGB masses from m 2 PGB /A 2 i = 0.04 — 0.2, so our S%PT expression 
must apply throughout this range. Generic discretization errors are of the size a 2 Aq CD , 
which is approximately 0.04 for 1/a ~ 2 GeV and Aqcd ~ 400 MeV, so they are compara- 
ble to the minimal m 2 PGB /A'^ and should be included at the same order. 2 Taste-breaking 
discretization errors, on the other hand, are caused by exchange of gluons with momentum 
it /a, and therefore receive an additional factor of ay(iT/a), so their size must be considered 
separately. 3 At the lightest quark masses, the lattice Goldstone pion mass is comparable to 
the mass splittings among the other PGB tastes: m 2 PGB / A 2 ~ a 2 al(q* = tt/o)A 2 ~ 0.04, 
where A is a QCD scale (distinct from Aqcd) which turns out to be around 1200 MeV. Thus 
they are not suppressed relative to pure 0(a 2 ) discretization effects by the additional powers 
of a, as naive power-counting would suggest, because of the large scale A associated with 
the taste-breaking process at the quark level jig . Standard S%PT only includes discretiza- 
tion effects - we now consider additional errors from perturbative operator matching, which 
depend upon ay(q*)- Generically, the choice of q* is process-dependent, and the value of 
ay(q*) ranges from ~ 0.3 — 0.55 for q* = n/a — 1/a at a = .125fm Thus both a/Air 

and a 2 must be included at lowest-order in our power counting. 4 

In light of this discussion, we adopt the following extended S%PT power-counting scheme: 



p 2 ~ m ~ a 2 ~ a 2 ~ a/An ~ a 2 , (15) 
where a 2 = ay(7i/a)a 2 . We account for the fact that a 2 a 2 terms in the action are enhanced 



2 Our generic discretization errors are not suppressed by a because we consider unimproved operators. This 
is in contrast to the mass and decay constant calculations of Ref . , which use improved operators and 

therefore have smaller generic discretization errors of around 2%. 

3 Taste violations are suppressed by a 2 rather than a because we assume an improved action. 

4 The quoted values of a indicate that it is perhaps necessary to include 0(a 3 ) errors as well; we assume 
that this is not the case. 
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numerically by including them at the same order as simple discretization effects. 5 In fact, 
while it may seem ad hoc, p 2 ~ m ~ a 2 a is the standard S^PT power counting scheme. It 
is simply not traditionally written as such because standard S%PT calculations have only 
included 0{a 2 a ) taste-breaking discretization errors from the action, and have therefore not 
needed to contrast them with pure 0(a 2 ) discretization effects. We also use conservative 
power-counting for the perturbative errors by assuming that 2-loop contributions are not 
significantly smaller than those from 1-loop diagrams. We emphasize that our scheme is 
phenomenologically based on the particular parameter values of current staggered simula- 
tions - simulations using significantly lighter quark masses or smaller lattice spacings would 
require a different scheme and result in a different S%PT expression for Bk ■ 

We must determine all of the contributions at NLO in our power counting to «Mi a t, 
including the effects of perturbative matching and discretization errors. We then find Bk by 
using A^iat in the ratio given in Eq. (|10|). In S%PT, Ai\ at is simply the matrix element of a 
sum of operators with undetermined coefficients: 

Mat = (k! p \C 1 x 1 x \K° 2p ), (16) 

so our goal is really to determine all operators in the chiral effective theory, O l x , that con- 
tribute at NLO to the above matrix element. Once we do so, the full NLO S%PT expression 
will contain the LO tree-level term, 1-loop corrections, and analytic terms of the same order 
as the 1-loop terms. 6 This task is not as daunting as it first seems given the extended power- 
counting scheme because it turns out that the only LO contribution to Bk is of 0(p 2 ). Thus 
operators of the following orders in our power-counting contribute to Bk at NLO: 

1-loop: 0(p 2 ), 0(m), 0(a 2 ), 0(a 2 a ), 0(a/An), 0(a 2 ) 
NLO analytic: 0{p 4 ), 0{p 2 m), 0{p 2 a 2 ), 0{p 2 a 2 a ), 0{p 2 a/4n), 0{p 2 a 2 ), 
0{m 2 ), 0{ma 2 ), 0{ma 2 a ), 0{ma/An), 0{ma 2 ), 
0(a A ), 0(a A a ), 0(a 2 /16vr 2 ), 0(a A ) . (17) 

In the rest of this section we give an overview of the operators that contribute to Bk at 

5 We assume that taste-conserving discretization errors are not enhanced since such behavior has not been 
observed in staggered lattice simulations. 

6 Note the distinction between Bk at NLO and Bk at 1-loop, which does not contain the analytic terms. 
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1-loop in S^PT. Appendix O contains the detailed operator enumeration as well as the 
determination of analytic terms. 

Clearly the lattice version of the continuum quark level operator Q s ^ a 99 ered ^ (jgj^ -will 
generate the dominant contribution to Bk , so we map Q s ^ a 99 ered on ^ c hiral operators first. 
This is relatively straightforward; we follow the graded group-theory method of Refs. 1qHi7| 
to insure that we include all possible linearly-independent operators. This results in two 
operators, Eqs. (jC8|) and (jC7|) . the first of which is simply the staggered, partially quenched 
analog of the continuum (27^, 1r) chiral operator, and the second of which arises because 
we use two sets of valence quarks. Because both chiral operators come from the same quark- 
level operator, their coefficients are related, and turn out to be equal. This is crucial because 
Bk in continuum %PT only depends on a single parameter Q], and our result must match 
the continuum one when a — > 0. Only the "standard" chiral operator contributes to Bk at 
tree-level, and its contribution is of 0(p 2 ). 

Next we determine the chiral operators which contribute to Bk because of perturbative 
operator mixing. Four-fermion operators which mix with Q s ^ a 99 ered are nQ ^ f as f e p are 
unaccounted for in Eq. (|1HJ) . so they introduce errors into Bk that are of 0(a/4%), as well 
as errors of 0(a 2 ). 7 We therefore map them onto chiral operators with two undetermined 
coefficients of 0(a/4%) and 0(a 2 ), respectively. The resulting eight chiral operators are 
given in Eqs. (jC16|) - (jC23|) . Taste P four-fermion operators are accounted for at 0(a/An) 
using the matching coefficients, but they still introduce errors into Bk at 0(a 2 ), through a 
new 0(a 2 ) operator given in Eq. (|C14J) . 

To match the continuum regularized operator, Q s ^ a 99^ed,cont ^ wg a j gQ neec i to include 
operators of higher dimension. From now on we will generically refer to this matching of 
continuum onto lattice operators as "mixing" . In particular, we must include all dimension 8 
operators that are in the same representation of the symmetry group that maps a hypercube 
onto itself. Because dimension 8 operators are explicitly suppressed by a factor of a 2 relative 
to dimension 6 operators, they can map onto chiral operators with coefficients of 0(a 2 ) that 
contribute to Bk at NLO. There turn out to be nine such operators - the same as those 
which arise from perturbative operator mixing. 



7 Recall that a /Air ~ a 2 in our power-counting because we make no assumption about the numerical size 
of the 0(a 2 ) matching coefficients. 
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Finally, four-fermion operators from the staggered action that arise due to discretization 
effects can enter time-ordered products with the Bk operator Q s ^ a 99 ered ? effectively correcting 
the Bk four-fermion vertex at 0{a 2 ). The diagram for this process is shown in Figure 01 
We must therefore determine the chiral operators which result from a combination of the 
Bk four-fermion operator and the inserted four-fermion operator from the action. We do so 
using the spurion method of Ref. Q|. This results in four new chiral operators, Eqs. (jC26|) - 
flC29|) . as well additional copies of seven operators which had appeared previously from 
operator mixing. 

In total, we find fifteen chiral operators; we will determine their contributions to Bk at 
1-loop in the following section. We note, however, that many of these chiral operators arise 
in more than one way and correspond to more than one quark level operator, so they ac- 
tually have more than one undetermined coefficient. We list the (9(a/47r), 0(a 2 ), 0(a 2 ), 
and 0(a^) operators and their coefficients in Table IIV1 For example, the operator 0^ A , 
which comes both from 1-loop perturbative matching and from insertions of operators in 
the staggered action has three coefficients, one multiplied by a/4ir, one multiplied by a 2 , 
and one multiplied by a 2 . The three coefficients can be separated in principle by carrying 
out simulations at more than one lattice spacing. We emphasize that they cannot be lumped 
together in a fit to multiple lattice spacings. In following sections, including the final ex- 
pressions for Bk at NLO, we only associate a single coefficient, C^,, which each operator. 
This is simply to reduce the size of the final expressions. One must combine the expressions 
with Table IIVI in order to make them complete. 

IV. B K AT 1-LOOP FOR 4+4+4 DYNAMICAL FLAVORS 

We first calculate Bk at tree-level, which receives a contribution from a single diagram, 
Figure EJa), and from a single chiral operator: 

°x = Yl (str(S9 M S t F 1 )Str(S^S t F 2 ) + p.c.) . (18) 

In the diagram, the Bk chiral operator is shown as two disconnected squares, each of which 
produces a kaon. This reflects the structure of Off, which has two supertraces in which the 
matrices F\ and F2 have the appropriate taste (£5) and flavor (sldl and s2d2) structures 
to produce a and a K%. is just the staggered analog of the continuum (27 l, 1r) 
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Bk operator - it reduces to the standard form after removing the taste structure (£5 — > £/) 
and the partial quenching (dl, d2 — > d and Str — > Tr). The tree-level Bk matrix element is 
simple: 

2C K 

M k° = — IT ^(/^IStr^TrFOStr^vrFs)!^ ) 



1 



/' 

32C K 



P 



32C* 2 



f 



2 x Vp' 



(19) 



where rn? xyp = fi(m x + m y ) is the tree-level kaon mass. The kaon B-parameter itself is just 
the ratio M.kI M-n^c-, where .M vac is given in Eq. Q. 8 Because the matrix element -M V ac is 
proportional to the square of the axial current, it is proportional to f 2 p 2 at tree-level: 

= \fml yp . (20) 

Therefore the kaon mass drops out of the ratio of matrix elements: 



LO 

= B , (21) 



and is just a constant at lowest order in S^PT. 

We are now ready to calculate Bk at 1-loop in S^PT. It is useful to simplify the form 
of Bk as much as possible diagrammatically before explicitly showing any expressions. In 
general, Mk receives 1-loop contributions from the five diagrams shown in Figures H^b)- 
(e) 3|. However, one can easily show that Figure ^f) cannot contribute in the staggered 
case because it is impossible to draw a corresponding quark-line diagram with two different 
external kaons. Moreover, because Bk is defined by the ratio M.k I ' M-vac, cancelations occur 
between matrix elements in the same manner as the kaon mass dropped out at tree-level. It 
turns out that diagrams W[b) and |TJc) cancel entirely in the Bk ratio, as we now show. 

The 1-loop matrix element M.k can be expressed as 

M K = \B Q fml yp {\ + X[Figs. UpbHc)]} + X'[Figs. md)-(e)] , (22) 



Because we work only in the staggered theory for the rest of the paper, we drop the superscript "staggered" 
from both operators and matrix elements in order to reduce clutter in expressions. 
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a. xi 

(a) (b) (c) 

XX XI u_ 

(d) (e) (f) 

FIG. 1: Tree-level and 1-loop contributions to Mk- In diagrams (a) - (e), one external PGB is a 
and the other is a K®. Diagram (b) generates kaon wavefunction renormalization. The circle 
represents a vertex from the LO staggered chiral Lagrangian. The squares represent an insertion of 
the Bk operator. Each box "changes" the quark flavor from d «-> s. Diagram (f) does not actually 
contribute in the staggered case because it has the wrong flavor structure to contract with two 
different external kaons. 

where X and X' denote the results of the specified diagrams and rn 2 xyp is the 1-loop kaon 
mass-squared. 9 The NLO expression for M. vac has the same form as the tree-level one: 

-^vac = \< y jly P , (23) 

except that both rr? and f xyp become the 1-loop quantities. Now consider diagrams [TJb) 
and^c). Visually, it is easy to see that they factorize - if one draws a vertical line between 
the boxes, the left halves of the two diagrams renormalize fx at 1-loop while the right halves 
just produce / at tree-level. It therefore seems natural that X[Figs. IH[b)- ( c )] should be 
related to the 1-loop renormalization of the kaon decay constant, and, if one works out the 
details, one finds that 

X[Figs. (c)] = 2^p° , (24) 

where the factor of two results from the fact that loop can be on either leg. Thus diagrams 

9 Note that X' is not within the curly braces because it comes from operators other than and is therefore 
not proportional to Bq. 
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(b) and (c) are exactly those necessary to turn the tree-level f 2 into the 1-loop /, 



2 .10 

xy P - 



M K = -B j% p m 2 xyp + X'[Figs. md)-(e)] , (25) 



and Bk only depends on diagrams [TJd) and[T^e): 

R i-ioop p , 3 ^[Figs.[H:d)-l(e)] 

B k = 5 o+8 fT—2 • ( 26 ) 

° Jxy P " L xyp 

The diagrams in Figs. E^d) and [IJe) receive contributions from all of the enumerated 



10 Note that the NLO analytic contributions to f xyp simply change the coefficients of the NLO analytic 
corrections to Bk- 
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Bk chiral operators, which we list here for convenience: 





= V (strfES^E'FOStrlSapStFj) +Str(S t a p EF 1 )Str(E , 3 ( ,EF 2 )) 

\ / 

A* 


gN 

X 


= V fstr(S9„S t F 1 S9 M S t F 2 ) + Six(E%EF 1 tf dj:F 2 )) 


1P 

X 


A* 

= Str(F 1 SF 2 S t ) 


of 


= Str^SFiS^sS^S^+p.c. 


of 


= Str(^EF 1 Et)Str(e 5 EF 2 E t )+p.c 


°? 


= Str(F 1/ SF 2/ E t ) 


of 


= ^ Str(F 1T EF 2T E t ) 


of 


= ^Str(^SF 1 St^SF 2 St)+p.c. 


of 


= ^Str^EFxE^Str^EF^) + p.c. 


of 


= ^ (str(F iy EF 2y E) + Str(F iy E t F 2y E t )) 
a* 


Of 


= ^ (str(F iy E)Str(F 2y E) + Str(F iy E t )Str(F 2V E t )) 

A 1 


of 


= ^Str(F iy E)Str(F 2y E t ) 


of 


= (str(F 1A EF 2A E) + Str(F lj4 EtF 2A Et)) 


of 


= (str(F 1A E)Str(F 2A E) + Str(F 1A Et)Str(F 2A Et)) 

A* 


of 


= ^Str(F 1A E)Str(F 2A Et). 



A' 



Appendix [0 explains the details of our operator notation, but we note here the important 
features. In general, each operator contains an Fi and an F 2 , which are taste matrices with 
different, nontrivial flavor structures such that produces the desired kaons at tree-level. 
In particular, F\ and F 2 with no additional labels are taste £5, whereas the remaining F's 
have "incorrect" tastes, e.g. Fit is taste The fact that many operators with tastes 
other than P contribute to Bk illustrates the importance of incorporating taste-symmetry 
breaking when calculating quantities with staggered fermions on the lattice. 
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(a) (b) (c) 

FIG. 2: Quark line diagram contributions to Bk at 1-loop. One external meson is a K® and the 
other is a K®. The two boxes (hexagons) represent an insertion of the Bk operator. Each box 
"changes" the valence quark flavor from (dl <-> si) or (d2 «-> s2). Each hexagon "changes" the 
valence quark flavor from (dl «-> s2) or (d2 «-» si). Diagrams (a) and (b) contribute to Figure [IJd) 
while diagram (c) contributes to Figure ^e). 

Figure C^d), m which one of the many Bk operators produces a four- meson vertex and 
two of the legs are then contracted to form a loop, corresponds to two quark-level diagrams, 
shown in Figures|2fa) and (b). The two boxes in each diagram are similar, but not identical: 
one changes a dl quark into an si quark while the other changes a d2 into an s2. Diagram 
(a) must be disconnected at the quark level because the two external kaons contain different 
kinds of valence quarks, so it only receives contributions from double-supertrace operators: 11 

D*M+ Di{q) -2D±(q) 



16(2Cf 4 


-Cf) 


f d'q 


/ 4 

16(2^4 




1 (2tt) 4 

r A 


f 4 




1 (2tt) 4 



\D v xx {q) + D V yy {q) -2D v xy {q)}. (28) 



The symbols D 1 , D v , and D A represent the taste /, V, and A disconnected propagators, 
respectively - we will show their explicit forms later as needed. Diagram Efb), on the 
other hand, is connected and only receives contributions from the seven single-supertrace 



11 Double supertrace operators 0^ p and cannot contribute to Bk at NLO because the loop mesons 
must be connected through a hairpin propagator, and can only be tastes /, V, or A. For a review of S%PT 
and brief discussion of hairpin (quark disconnected) propagators see Appendix IbI 
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operators 12 : 



nK fB 1 



B' 



d 4 q (2q 



plB n BB' 



B,B 



£ 

B,B' 



' f 4 



d A q ri 

(2^\ 3 



(27T) 



- 


1 




+ m 2 , 

X VB' 




1 




X VB' 




1 


g 2 


+ m 2 

X VB' 



(q 2 + p 2 ) 



1 



1 



«„/ + ^ + m TO S ' 



1 1 



4Lg 2 + m 2 XR; g 2 + ^ ys; 



(29) 



where the coefficients f B ', g BB \ and h BB ' are matrices in which B = I, P, V, A, T labels the 
contributing operator and B' = I, P, V, A, T indicates the loop meson taste: 



/ 



B' 



8 [114 4 6 



(30) 



BB 



I 


1 


1 


-4 


-4 


6 \ 




1 


1 


4 


4 


6 




-8 


8 


-16 


16 







-8 


8 


16 


-16 





\ 


12 


12 








-24 / 



(31) 



,BB' 



128 



(32) 



/o \ 
110 


\0 6 6 12 / 

The sum over B' accounts for how many times a particular loop meson taste contributes to 
the above expression with a given coefficient - this is important because a staggered meson's 
mass depends on its taste as well as its flavor. 13 The pattern {1, 1, 4, 4, 6} in Eq. (f3T)|) counts 
the number of lattice PGBs with each taste - thus the coefficient f B ' averages over all PGB 
tastes. Note that the rows of zeros in Eq. (|32|) have no deep physical meaning, but simply 



correspond to the absence of certain operators, e.g. h 0B = (0, 0, 0, 0, 0) indicates that O 21 



12 Recall thai Of and Of have the same coefficient, which we call Cr; . 



13 



Sec Ref. |£( for the staggered meson masses at tree-level. 
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does not exist. Nevertheless we include them so as to allow both B and B' to run over the 
same set of indices. 

Figure ^e), in which an operator from the staggered chiral Lagrangian generates the 
four-meson vertex and one of the Bk operators is inserted in the loop, corresponds to the 
third quark level diagram, Figure|2fc). Because the LO staggered chiral Lagrangian contains 
three types of operators - the continuum two-derivative term, the continuum mass term, and 
the 0(a 2 ) potential - there are many contributions to this diagram. Nevertheless, although 
the staggered potential contains both single- and double-supertrace operators, we only show 
a connected vertex in Efc). 14 This is because it turns out that, for all disconnected 4-PGB 
vertices from the LO staggered potential in which two of the legs have the desired flavors 
of external kaons, the various Wick contractions cancel and the vertex is identically zero. 
Thus we need only consider connected 4-PGB vertices. Because the two external kaons 
have entirely different flavors, they cannot share any quark lines without a Bk operator 
insertion (box or hexagon), so Figure Efc) shows the only possible quark flow. However, 
it is now easy to see that the Bk operator must change both the quark flavor and the 
valence set, so we use hexagons rather than boxes to indicate this. Finally, at the 2-PGB 
level, double-supertrace Bk operators (by construction) change flavors within a valence set, 
whereas single-supertrace operators change flavors between sets of valence quarks. Thus the 
last diagram only receives contributions from single-supertrace operators: 

7l, 



^ = E/ (^[-cf/ B '« 2 +E(4V fl '-cf^' 

1 



X 



r p 2 + q 2 - 2m 2 xyp + m 2 xyBi 



6/ 4 



q 2 + m 2 



(33) 



xy B ' ■ 

Note that it depends on the same coefficients, /, g, and h as - this will allow us to 
combine the two in the following section. 

V. NEXT-TO-LEADING ORDER B K RESULTS 

The 1-loop matrix elements determined in the previous section apply to a PQ theory 
with 4 tastes of each sea quark flavor. To go from four tastes per sea quark to one taste per 



The staggered potential is shown explicitly in Appendix [5] 
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sea quark (from the so-called "4+4+4 theory" to the "1+1 + 1 theory" when N sea = 3), we 
must multiply every sea quark loop by 1/4. At first glance, the quark- level diagrams which 
contribute to Bx at 1-loop, Figs. |2Ja)-(c), do not appear to contain any sea quark loops. 
However, sea quark loops implicitly enter the disconnected hairpin propagators when one 
rediagonalizes the mass matrices of the flavor-neutral, taste V, A, and J, sectors. In general, 
the forms of quenched, PQ, and full QCD matrix elements only differ for diagrams with sea 
quark loops, so the final expression for Eq. (|28[) will change in the three cases, but those for 
Eqs. (129(1 and ()33() will remain the same. 



A. B K at NLO in the PQ Theory 

We first consider Bk in the 1+1 + 1 PQ theory (m u ^ ma ^ m s ), as it is the most general. 
We label the valence quarks by x and y, and reserve the labels u, d, and s for the sea quarks. 
Note that the sea quark masses only show up in the PQ results implicitly through the masses 
of the flavor-neutral PGBs: 7r°, n, and rf. Thus the 1 + 1+1 result can be turned into the 
2+1 {m u = ma ^ m s ) result simply by changing the expressions for the 7r°, t], and r{ masses. 

The matrix elements that come from connected quark-level diagrams are most simple to 
calculate, so we consider them first. Neither Fig|2(b) nor Fig. Etc) contains sea quark loops, 
so we do not have to insert any factors of 1/4 by hand - we need only perform the integrals. 
To simplify the resulting expressions, we use a condensed notation for the chiral logarithms. 
As in Refs. g| and we define two functions, £ and I: 



d 4 q 1 1 



(2vr) 4 q 2 + m 2 ' 16tt 2 
d 4 q 1 1 - 



(2?r) 4 (q 2 + m 2 ) 2 ' 16tt 2 
The remaining integrals can all be expressed in terms of I and I: 



m 2 ) (34) 



m 2 ). (35) 



" m V 2 ) (36) 



(2tt) 4 q 2 + m 2 16tt 2 
d 4 q q 2 1 



m 2 ) - m 2 £{m 2 ) (37) 



(2tt) 4 (q 2 + m 2 ) 2 16tt 2 

' m 2 £{m 2 ) - 2£(m 2 ) ) . (38) 



d q q 4 

(2vr) 4 {q 2 + m 2 ) 2 = 16^ 
Because finite-volume (FV) corrections to Bk only alter the chiral logarithms by turning 
integrals into sums, one can use either the infinite volume or FV expressions for £ and £ as 
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desired: 



Om 2 ) = m 2 (\n^— + <5f y (mL)), 5? v (mL) = -i- ^ 

,2 



\r\ 



e(m 2 



(ln-J- + l) +5 3 Fl/ (mL), 5 3 Fy (mL) = 2^ifo(|r1mL), 



(39) 
(40) 



where L 3 is the lattice spatial volume and K and K\ are Bessel functions of imaginary 
argument. 15 Thus our final expression for Bk will be completely general. In terms of I 
and £, the 1-loop matrix element contribution to Bk from connected quark-level diagrams, 
M®+M®, is 



_J X^)c K f B ' 

64tt 2 / 4 \ x J 



m \pA m \j + m YpA m Y n ) - 2m ly B A m ly n ,) 



xy B > 



[C 1B g BB ' + C 2 x B h BB ' 



£(m 2 x J+£(m^,)-2£(m 2 xyB/ ) 



i{m\J + i{m\J + 2£(m 2 xyB/ ) - 2m 2 xy J{m 2 xy J 



xy B r\ xy B , 



- £ ( c f s BB ' - c i" hBB ' 



xy P 



2i(m 



2 

xy B , , 



We use a condensed notation for the flavor-neutral mesons, in which m 2 x 



(41) 



mL and so 



forth. Note that the first two lines vanish when m. 



m y . Note also that the second and 



third terms in M. con n are proportional the mass of the external kaons because we have let 



2 2 

p = —m 



xy P - 



Next we consider the matrix element that comes from the disconnected diagram in 
Fig. |2Ja), Eq. (|28|). which contains tastes /, V, and A hairpin propagators. Because the 2- 
PGB disconnected vertices are dimension 2, an additional hairpin vertex can be canceled by 
an additional propagator in the S^PT power-counting, so one must resum the flavor-neutral 
propagators with all possible numbers of sea quark loops between the external hairpins. 
The staggered, flavor-neutral, full propagators were determined in Ref. j^] using the general 
method outlined in Appendix A of Ref. j^j, but adding factors of 1/4 for every sea quark 
loop in the series. We will just quote the results. 

Generically, hairpin propagators contain poles at the mass eigenstates (it , 77, and 7/) and 



is 



Note that we use dimensional regularization and the same rcnormalization scheme as Rcfs. 



HQ. 
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the flavor eigenstates (X and Y). 16 These multiple-poles can be rewritten as sums of single 
poles times their residues, e.g.: 



XV ' V (q 2 + m 2 Xv )(q 2 + m^ v ) (q 2 + m^q 2 + m 2 v )(q 2 + m 2 ^) 1 ) 



_ fl2 1 (g 2 + m£ v )(<Z 2 + m2, v )(g 2 + m 2 ^ 

= -a% ^ 2 1 2 Rff\{m x ,m Y ,m n o,m v ,m vl }]{m u ,m D ,m s }),(43) 

j=X,Y^°,ri,v' jv 

where Rff is the residue of the single pole at q 2 = —m 2 v , e.g.: 

^[5,3] = Ky - m ly)« ~ ~ m Xy) (44) 

X ^ (mfr y - ) (m 2 „ - m\ y ) (m^ - m\ y ) (m 2 ^ - m 2 ^ ) ' 

8| and 0, which is most clearly described in the Appendix of 
0|. The expression in Eq. (|42|) is valid for 3 sea quarks with any combination of masses; an 
analogous expression holds for D^ PQ (q). Because the overall flavor-taste singlet decouples 
from the PQ theory, when m — > oo the expression for D^ >( ^(q) simplifies further: 



We use the notation of Refs. 



4m 2 1 (g 2 + m 2 Ui ) (g 2 + m 2 Dl ) (g 2 + m% ) 

3 



D I,PQ (q) = - K - UlJKH DlJK ^ blJ (45) 

(g 2 + m 2 Xi )(q 2 + m 2 Yi ) (g 2 + m 2 ? )(g 2 + m\){q 2 + m 2 , ) 1 } 



4 1 (q 2 + ml^q 2 + m 2 Di )(q 2 + m 2 Si ) 



(46) 



oo 3(g 2 + m| / )(g 2 + m 2 7 ) (g 2 + m 2 o)(g 2 + m 2 n ) 

-\ £ ^- L ^<' 3] (W,m y ,m 7r o,m,};{m c/ ,m D ,m 5 }).(47) 



Some of the taste / disconnected propagators are multiplied by g 2 , so the residues of their 
poles are changed: 

q 2D l y PQ (q) ~ \ Y) 2 I 31 2 Rff\{rnx,m Y ,m T o,m ri };{m u ,m D ,m s }), (48) 
3 z — ' q z + rri: J1 

j=X,Y,7r°,ri * 31 

where "~" indicates that this relationship only holds within the integral. 
Using the definitions of the disconnected propagators it can be shown that 

D xx + D yy - 2D xy = (m 2 x - m Y ) 2 -^-^{D xy }. ( 49 ) 



16 The staggered factors of 1/4 shift the mass eigenvalues but do not introduce additional poles. Expressions 
for the 7T°, 77, and rj' masses in a theory with 2+1 sea quarks are given in Ref. The generalization to 
the 1+1+1 case is straightforward. 
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This greatly simplifies the 1-loop matrix element contribution to Bk from disconnected 
quark-level diagrams: 

■>K 



(2Cf + Cf)a 2 5' 2 2 2 9 9 / ^ 2, p [ 5l 3K\ 



j=X,Y,Tr°,r),v' 



(2Cf + Cf )a% 2 2^99/^.. MpM\ r ^ 



vr 2 / 4 

(2Cf + Cf)a 2 ^^ 2 ^ 2 , 2 9 9 

y v v j=X,Y,ir ,r),r]' 

where we have set p 2 = —m 2 yp and the arguments of the i?'s are as in Eqs ()43)1 and (j47jh This 
compact notation emphasizes that the disconnected matrix element vanishes quadratically 
as {m x — m y ) — > 0. 

Finally, at NLO in a PQ theory with 1+1 + 1 sea quarks is 



'X "^o I f2 ™2 /2 r f2 r ° /2 ™2 r ^ /2 



, 3 M conn + M d £ c A m xyp (rn x -m y y (m u + m d + m s 



xyp xyp J xyp J xyp J xyp xyp J xy P 

(51) 

where £>o is given in Eq. (j2"Tj) and .M conn and M.j® c are given in Eqs. (|4~Tj) and (jBTljl . Here 
we have reintroduced the analytic terms from Appendix IC 31 We have checked explicitly 
that all renormalization scale dependence in the logarithms can be absorbed by the analytic 
terms. 



B. B K at NLO in the Quenched Theory 

We now find Bk in the quenched theory. Because of the absence of sea quark loops, 
additional factors of 1/4 to reduce the number of tastes per flavor are unnecessary. Moreover, 
the flavor- neutral propagators need not be resummed, so their forms are quite simple, e.g.: 

D VQuench^ = - a 2 5 > 1 (52) 

y (<? 2 + m 2 Xv ){q 2 + m\ v ) 

and likewise for D x ^ uench (q). However, the flavor-taste singlet (tj'j) does not decouple in the 
quenched theory, so we cannot take mo — > oo: 

n i, Quench (n] 4 m 2 , + aq 2 

xy {q) ~ ?>{q 2 + m 2 Xi ){q> + rn Y y 1 ' 
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where a is an additional chiral parameter that is only present in the quenched theory, not 
the strong coupling constant. As before, the nondegenerate propagators can be further 
simplified into sums of poles times residues: 

1 1 



D VQuench^ = _ Q 2£/ 



Y v Ay 



q 2 + m Xy q 2 + m\ v 



-4 \ mi — amis m 2 , — ami 



D I, Quench, , = I Zl ) t "*° ZZZl _ '"° _ "'"^ (ZA) 

V ' V 3 )m 2 Yi -m 2 Xi l(q 2 + m 2 Xi ) {q 2 + m 2 Yl )\' K } 

Because the connected quark-level diagrams that contribute to Bk at 1-loop, Figs. I2jb) 
and (c), do not contain sea quark loops, the corresponding matrix elements remain un- 
changed from the PQ theory. The matrix elements from disconnected quark-level diagrams, 
after setting p 2 = —m 2 xyp , are 

^ST* = ^«-O a a^a^{ E/K)(<, + -l)K-^)< ] } 

(2Cf + C^)a 2 5' A 2 2 9 9 / ~ [2i0] 

(2Cf + Cl A )a 2 5' v 2 2 9 5 / ~ [2>0] 1 

+ ^3 K " ™>v) dmY^i 22 ^ m iv) R\ v Y (55) 



where 



R l x7 = R l x? m r}) = 7Z-2 = "^v ■ (56) 

Therefore B^- at NLO in the Quenched theory is 



\m Yv -m Xvi 



R Quench R , 3 ( M conn + Mf s e c nch A m 2 xyp (m x - my ) 2 \ 

\ Jxy P ""xyp Jxyp J xy P Jxy P " u xyp J 

where B is given in Eq. (J2T|) and A4 conn is given in Eq. (jUJ). We emphasize that the 
undetermined coefficients, C l and A-C, are different in the quenched theory than in QCD, 
and that there is no D term. 



C. B K at NLO in the Full (2+1) Theory 

Finally we determine Bx in full QCD, i.e. m x = m^ and m y = m s , with 2+1 sea quarks 
(m u = m^ m s ). Although one can, in principle, just take the limit of the final PQ 
expression, Eq. (|51j). it is easier to start from the expressions for the 1-loop matrix elements 
in terms of integrals. 
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Because QCD is a physical theory, the disconnected propagators only have single poles. 
It is easy to show this explicitly for the 2+1 theory, in which m 2 ^ = mfj = m 2 D and parts 
of the numerator and denominator cancel. For example, the taste V (and A) disconnected 
propagators are: 

1 



D 



D 



v dd 2+1 (q) 

V, 2+1 , 



(?) 




{m 2 Sv -m 2 Dv ] 



_q 2 + m 2 Dv \ {m 2 v - m 2 Dy ) (m^ - m 2 Dv ) 



+ (D -> 7] -> rj[ -> £>)(58) 



V, 2+1 



(?) 



2 r/ 



(77 «-> 77') 



(59) 
(60) 



where the arrows indicate all possible permutations of masses. As in the PQ theory, the 
taste singlet 7/ decouples, so we can use the explicit large-mo expressions for the masses, 



m u , = m 



m 



m 
m 



Dj 



Di 

2m\ 



m 



o- 



(61) 

(62) 
(63) 



to drastically simplify the taste / disconnected propagators: 



D 



1,2+1 

dd 



D 



7,2+1 



D 



1, 2+1 

ds 



(?) 



(?) 

(?) 



mo^oo q 2 + m 2 Di 3 q 2 + m 2 m 
1 8 1 



-4- 



q 2 + m 2 Si 3q 2 + m 
4 1 



2 
m 



(64) 
(65) 
(66) 



mo^oo 3 q 2 + 

In this section we only show the 1-loop matrix elements that come from disconnected 
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quark-level diagrams: 



M%1 



vr 2 / 4 

72V 



m 2 dsp + m 2 Di )i(m 2 Di ) + 2(m 2 dsp + m 2 Si )£(m 2 Si ) - 3(m 2 dsp + m 2 j£(m 

3V\„2r/ 



(2Cf + Cf)a 2 5 A 



£(m 2 

2 



Da ' 



\ m S A - m D A , 



Sa> 



(m 2 A -m 2 DA )(m 2 , A -m 2 DA ) 
( m D A ~ ™ 2 Sa ) 



(ml A -m 2 SA )(m 2 , A -m 2 SA ) 



(2C 2A + q A )a 2 5' v 

2 f4 



vr 2 / 



^Kv) (m 2 _ m 2 ) 



2 ' 



[m Sv -m Dvi 



(m 2 v - m 2 Dv )(m 2 v , v - m 2 Dv ) 

{m 2 Dv -m 2 Sv ) 
(m 2 vv -m| v )(m^ -m| v ) 



+ (?7 <-> r/' 
+ (D -> 7/ -> rj[ -> D) 



+ (5 — > 77 — > 77' — > 5) 



(67) 



The matrix elements which come from connected diagrams are minimally changed from the 
PQ ones in that x — > <i and 7/ — > s. Therefore £>x at NLO in full QCD with 2+1 sea quarks 
is 



3 

+ - 



.4 



B- 



m 



dsr +C 



(m d - m s ) 2 ^ (2m d + m s 



8 \fds P fdsp fds P m dsp fdi 

where B is given in Eq. (|2*T|) and M. conn and is given in Eq. (jUJ). 



(68) 



VI. GUIDE FOR LATTICE DETERMINATION OF B K 

Any lattice determination of Bk necessarily entails two steps: first, fit the expression 
to the lattice data and determine the unknown coefficients, and, second, use the resulting 
function to extrapolate Bk to its value at the physical quark masses in the continuum limit. 
The large number of undetermined constants in the staggered %PT expression for Bk at 
NLO could make the fitting procedure quite difficult, and certainly necessitates a lot of data. 
Our goal is therefore is to find alternative ways to access some of these coefficients. Once 
measured, they can be used in the full NLO S%PT Bk fit. 
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We begin by reviewing explicitly all of the parameters in the expression for Bk at NLO. 
We presume that Bk will be fit as a function of the kaon mass squared, so we do not consider 
the latter to be a fit parameter. In order to properly count the number of undetermined 
coefficients, we must combine the PQ expression for Bk, Eq. 1)51)). with Tables HVI and IV! 
Only thirteen of the fifteen operators listed in Eq. (J27)) in fact appear in Eq. ()51|). Operators 
0^ p and C^ T turn out not to contribute to B K at NLO. Operators O 2 ^ and (and 
similarly and 0^ ) turn out to enter the expression for Bk as a single linear combination, 
and can therefore be associated with a single coefficient. Thus we can see from Table IIVI 
that there are 23 coefficients associated with NLO chiral operators, in addition to . 
Similarly, Table |3 shows that there are 6 linearly-independent analytic terms. Finally, 
certain parameters appear in all S%PT calculations: the four tree-level mass splittings, the 
two hairpin parameters, and f X y P - Thus the expression for Bk at NLO is given in terms 
of 37 undetermined coefficients. Fortunately, some of these parameters have already been 
determined on the MILC 2+1 dynamical field configurations Q|. In particular, the taste 
Goldstone meson decay constant, as well as the mass splittings between the other PGB 
tastes, are already known. 17 However, this still leaves 32 undetermined fit parameters. 

An obvious way to reduce the number of parameters is to perform a fit to a single 
lattice spacing. This reduces the coefficients in Table IIVI to one per operator (or per linear 
combination), i.e. from 23 to 9, and the analytic terms in Table fVl from 6 to 4. Clearly 
one must eventually fit to multiple lattice spacings. We suggest, however, that is preferable 
to fit Bk in a multi-stage process. First perform separate NLO Bk fits at multiple lattice 
spacings with this reduced set of parameters. As we will show, a single lattice spacing 
fit of Bk, in combination with tree- level fits of other matrix elements, can be used to 
separate almost all of the perturbative and discretization errors from the desired continuum 
result. Second remove these errors at each lattice spacing. Third fit the results after error 
subtraction simultaneously to multiple lattice spacings using the appropriately simplified 
S%PT expression. Last send the remaining errors, which have now been obtained, to zero, 
and calculate a value for the continuum parameter Bk- 18 

17 We note that the vector hairpin is poorly determined from the MILC data, and that effects of the vector 
and axial vector hairpins are, in practice, difficult to disentangle in S%PT fits. 

18 The explicit expression for Bk in terms of Bk is given in Eq. (|82f> . 
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We now discuss step one: fitting the NLO expression for Bk at a single lattice spacing. 
Because this still involves 16 parameters, it is qualitatively useful to first consider the case 
of degenerate valence quarks, i.e. m x = m y . Recall from the previous sections that the NLO 
expression for Bk greatly simplifies in this limit. The matrix element from the disconnected 
diagram vanishes, thereby removing coefficients - C^ A and the hairpin parameters, as do 
half of the terms in the matrix element from the connected diagrams. The degenerate mass 
case is therefore the appropriate testing-ground for whether the NLO expression describes 
the data at all, i.e. if one is in the chiral regime, before moving on to the more complex, but 
physical, situation of m x ^ m y . 

In fact, the degenerate mass case is also quantitatively useful. Recall that the domi- 
nant contribution to Bk at NLO comes from the tree-level matrix element of , and is 
proportional to C^: 

/JFo i . i'+oc+o v>\ l t^o \ c 

6-£ + .... (69) 



(K% | Ox-^taste P) | K%) = Cj 
-M vac / 4 



Because we wish to determine Bk (and consequently C£) to NLO accuracy, the above ex- 
pression is insufficient. 19 An NLO fit to Bk with degenerate valence quark masses, however, 
can potentially be used to determine to NLO accuracy. In this limit, Bk at NLO has 
the following form: 




flB f>2B 



J x vp b' b * 



+ A W~ W~ sp ' m 

J xyp J xyp J xyp 



It depends on C£ , now at NLO, three analytic terms, and five linear combinations of other 
operator coefficients, ^ B {C^g BBl - C^ B h BB ') for B' = I,P,V,A,T. As we show below, 
these five combinations can be determined to the required accuracy using other matrix 
elements. Once they are found, we would like to use a fit of Eq. (|7U|) to extract to NLO 
accuracy. Unfortunately, this turns out not to be possible; separation of the NLO part of 



19 In fact, at first glance, we also need / at NLO. However, this is really unnecessary because enters 
both the LO and NLO contributions to Bk as the combination C^// 4 . Thus we can treat the quantity 
(CyV/ 4 ) as the undetermined fit parameter, and likewise for the other NLO coefficients. 
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C£ from the analytic term A requires a simultaneous fit to multiple lattice spacings. We 
can, however, consider their sum, (C£ / f A + 3v4/16/J yp ), to be a single fit parameter that 
can be determined NLO accuracy. We can also determine the remaining analytic terms, B 
and D, using Eq. (J70)l at a single lattice spacing. 

The other simple handle that one has in lattice simulations, besides the quark masses, is 
the flavors and tastes of the external states. At tree-level, double supertrace operators have 
nonvanishing matrix elements between a K® (sldl) and a (s2d2). In contrast, single 
supertrace operators have nonvanishing matrix elements between "mixed" kaons: sld2 and 
s2dl. Both single and double supertrace operators can have nonvanishing matrix elements 
between all five possible tastes: /, P, V, A, and T. Thus, by evaluating the matrix element of 
C*i-2oop(taste P) between different external states, we can naively hope to determine as many 
as nine independent coefficients for use in the larger NLO S^PT fit. In essence, while the 
two sets of valence quarks and multiple tastes introduce many new operators into staggered 
xPT, they also allow the independent determination of more operator coefficients. 

First consider the double supertrace operators, of which there are four: O^f , ^x*' 
and . Clearly, to extract the coefficient of (or ) we need to choose taste vector 
(or axial) kaons for our external states. The resulting matrix elements are: 

(KfylO^itasbe P)\K° Y ) = — (2Cf -Cf) (71) 

(^Id-^Ctaste P)\K° 2A ) = -^(2Cf - Cf ), (72) 

so we can only determine the linear combination of coefficients (2C^ — ). However, in 
the 1-loop contributions to B^, these coefficients enter in a different linear combination: 
{2C^ + C^). Thus we cannot in any simple way separately extract the coefficients of the 
two supertrace operators. Nevertheless, calculation of these matrix elements should still be 
done, as it provides an estimate of the size of the NLO coefficients. Moreover, because the 
matrix element of 0i_£ oop (taste P) between taste / and T kaons should be zero at this order: 

(tffjlOi-to^taste P)|if 2 ° 7 > = (+NNLO) (73) 
(^lOi-too^taste P)\K° 2T ) = (+NNLO) , (74) 

they can be used to gauge the size of the higher-order terms, thereby testing the validity of 
our phenomenological power-counting scheme and determining whether or not it is necessary, 
as is done in Ref. 0, to introduce a small subset of NNLO terms. 
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Next consider the single supertrace operators. Extracting many of their coefficients turns 
out to be feasible, although not completely trivial. There are eight such operators - Of 
Of, Of, Of, Of, Of, Of, and Of - and five possible matrix elements. We can 
combine them into a single equation using the same coefficients as in Eq. (J7U|) : 



(iv|^r|Oi-/o Op (tasteP)|K 2 ™ x ) 3 



nK fB'2 



B ^ 



(75) 



where 



N B ' = ( -1 -1 4 4 -6 ) • (76) 



The coefficient can be separated through its quark mass dependence, but this is not useful 

A. 

because it is only determined at LO in this tree-level matrix element. However, the numerical 
values of the five linear combinations of coefficients corresponding to B' = I, P, V, A, T are 
needed, along with the degenerate mass Bk result, Eq. (J7U|) . to determine the combination 
of coefficients (12C^ / f 4 + 3A/8f 2 yp ) to NLO accuracy. The five matrix elements in Eq. (f?3|) 
also aid in the true, m x ^ m y , E>k determination because we can solve for five coefficients 
in terms of the remaining two, and substitute the resulting functions into the NLO S^PT 
expression before fitting. 

Thus, after exploiting both degenerate mass and mixed-flavor matrix elements, we are left 
with an NLO fit at a single lattice spacing to only seven parameters: 5' v , 5' A , (2Cf + Cf ), 
(2Cf + Cf ), C, and two additional NLO coefficients of our choice from Eq. ()75|). This 
should certainly be feasible. Once this is done, the next step is to remove as much of the 
discretization and perturbative error as possible before combining the data from multiple 
lattice spacings. We do this by setting all coefficients that are multiplied by a 2 , a 2 , a/An, 
and a 2 to zero. 20 One can see now why it is preferable to fit to one lattice spacing at a 
time - it avoids having to give each operator in Table I1VI a separate coefficient for each 
expansion parameter. Doing so would clearly be a waste of effort, since one only needs the 
total error, not the individual coefficients, in order to remove it. We emphasize, however, 



20 This is obviously the correct way to remove discretization errors. The right way to remove perturbative 
errors, however, is less apparent, since we have only partially matched the continuum operator at 1-loop. 
We will later show that it is, in fact, correct to simply set their coefficients to zero. 
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that there is one coefficient that cannot be removed at this stage - A. It is tied up in the 
linear combination (12C^// 4 + 3A/8f 2 yp ). We must therefore remove all of the coefficients 
of operators in Table HVI and then perform a new fit to multiple lattice spacings. 21 Such a 
fit will only contain the continuum parameter and A, the latter of which must now be 

A 

separated into three pieces of 0(a 2 ), 0{a 2 a ), and 0(a 2 ). It therefore has four parameters, 
which is a clear improvement over the 37 parameters that would have originally been needed 
in a fit to multiple lattice spacings. Once we have separated the three contributions to A 
from the coefficient at NLO, we will remove them as well. 

Finally, assuming that we have successfully carried out the fitting procedure outlined 
above, we are ready to turn the result into a continuum value for Bk- We therefore review 
the procedure for matching lattice and continuum operators. Along the way, we clarify some 
of the points and justify some of the assumptions made previously. 

After removal of discretization errors, multiple lattice regularized operators match onto 
a given continuum operator in the following manner: 

O cont (/i) = Z(fi,a)0 Ut (a), (77) 

where (9 cont (/i) is a single operator, O lat (a) is a vector of operators, and the vector Z relates 
the two. The 1-loop expression for a) is typically written in the following way: 

Z(n, a) = 1 + — c(/ia) , (78) 

47T 

where c is a vector that contains the 1-loop matching coefficients. However, Eq. (|78|) clearly 
has a problem: what is the appropriate scale at which to evaluate a? This ambiguity can 
be avoided by using the exact perturbative formula for Z given by Ji in Ref jlo| : 

Z( ft a)=exp^-y o ^^-^jr^-J^g^j, (79) 

where T g i indicates a (/-ordered product 22 and 7 and (3 are the anomalous dimension and 
beta function, respectively. 23 Equation (fTTIj) is easy to understand physically; it shows that 



21 Strictly speaking we should fit to the quantity Bk, which is independent of lattice spacing, as we explain 
below. 

22 Note that there is no T g n because we are considering the case in which there is only one continuum 
operator and no mixing on the continuum side. 

23 Note also that the coupling g{a) can be chosen to be a continuum-like coupling such as MS j2p| . 
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matching is a three-step process. First run the lattice operators calculated at a scale q* = 1/a 
up to a = 0. Next match the lattice operators at this scale onto the continuum operator at 
a scale \i = oo. Finally run the continuum operator down to a scale fi, which is typically 
2 GeV. This procedure cleanly separates the lattice scale, a, from the continuum scale, \x. 

Mathematically, Eq. (|79|) can be used to derive the following relationship between the 
continuum and lattice operators: 



a ( ? *)-70/2ft 

where 



1 + ^1 j cont + (a 2 ) O cont (/i) 
1 + ^ Jcont + O(a') 

4:11 



1 + ^flc(q*a) + 0(a> 
Air 



C lat (a) , (80) 



T - A7 ° 71 (81 ^ 



and the vector c contains the standard 1-loop matching coefficients. This formulation of 
perturbative matching is extremely useful. Because the matching coefficients only appear 
on the "lattice side" of the equation, and no longer depend on the continuum scale one 
can choose any reasonable scale at which to calculate the matching coefficients - a natural 
choice is q* = 1/a. One must, of course, then evaluate all a's in the NLO S%PT expression 
for Bk that arise due to truncation of perturbative matching factors at this same scale. 

Equation (|8T)|) also justifies our treatment of errors associated with perturbative matching. 
Recall that the standard matching procedure is to only match to taste P operators at 1-loop. 
In terms of the above matching expression, this means that one calculates the matrix element 
of [1 + c(a/4-7r)](9p t , where the vector O l p contains only taste P operators. However, this is 
not correct, as lattice operators of all tastes are required to correctly match the continuum 
operator. We account for this error by including "wrong-taste" operators in the S%PT fit 
with coefficients of 0(a/Air). Conversely, if we were to correctly match onto the continuum 
Bk operator at 1-loop, we would not need to include such operators. Thus we conclude that 
in order to extract the correct result for Bk in the continuum we must set the coefficients 
of 0(a/An) in Table IIVI to zero. One can use the same logic to show that all of the 
perturbative matching errors should be dealt with in this manner, i.e. we must also set 
the 0(a 2 ) coefficients in Tables IIVI and IS to zero. This result is fortunate, as it allows 
treatment of all errors - both from discretization and perturbative matching - in the same, 
simple manner. 



30 



Finally, Eq. (|80|) shows that, once one has successfully fit the S%PT expression for Bk and 
removed the discretization and perturbative matching errors at each lattice spacing, con- 
structing the renormalization group invariant quantity, Bk is trivial. This is because, to 



1-loop order 



2l| 



so 



B K = «(/i)" 7o/2/3 ° 



B K = a(q*)- Jo/2f3 ° 



a(/i) 

J- t . '-'cont 

Air 



B% S (»), (82) 



, ^cont 
47T 



7~>lat, subtracted /oo\ 
K ■> K 66 ) 



where £>^ t ' subtracted j s lattice value for Bk at the scale q* with discretization and pertur- 



bative matching errors subtracted. Conversely, we can use Eq. (|83j) to remove the remaining 
errors from Bk in the following way. Once we have values for Bk with all errors but those 
from A removed at each lattice spacing, we turn them into values for Bk by multiplying by 
a(g*)~ 7o//2/3 ° [1 + J corit a(g*)/47r]. We then simultaneously fit B K at multiple lattice spacings 
to separate the continuum parameter from A. Next we subtract A to get a continuum 
value for Bk- Last we extrapolate the continuum result to the physical quark masses using 
the continuum %PT expression. 

Even with the implementation of all of these suggestions, fitting the staggered Bk data 
will be highly nontrivial. One would clearly like to somehow decrease the number of op- 
erators, or at least undetermined coefficients, that contribute to Bk at NLO. One way to 
do this is by using improved links. This drastically reduces the 1-loop perturbative mixing 
between 0\_{ oov (taste P) and certain wrong-taste operators. In particular, for the case of 
hypercubic fat (HYP) links, the 0(a/Air) mixing with taste S and taste T op erators is so 
small that we can consider it to be of higher-order in our power-counting |l4j. This leaves 
the operators O^ 1 and (9* T each with only a single coefficient of 0(a 2 ), thereby reducing 
the total number of fit parameters by two. 24 

The biggest source of potential improvement, however, is in better perturbative matching. 
A good first step would therefore be to match to all tastes of lattice operators at 1-loop, since 
the requisite matching coefficients are known. While it has been thought that "wrong taste" 
operator contributions to Bk would be less important than those from taste P operators, our 



We presume that a similar simplification occurs for Asqtad links, although their perturbative matching 
coefficients have yet to be determined. 
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power-counting and operator enumeration show that this is not the case in general. Complete 
1-loop matching would eliminate the 6 coefficients of 0(a/47r) in the first column Table ITVl 
Moreover, 2-loop matching of all tastes, or better yet fully nonperturbative matching, would 
completely remove the first column in Table II VI and eliminate the operators O} 1 and O x J 

X, A. 

altogether at this order, reducing the total number of fit parameters from 37 to 24. 25 

VII. CONCLUSIONS 

We have calculated Bk to NLO in S%PT for quenched, PQ, and full QCD, the form of 
which is necessary for correct continuum and chiral extrapolation of staggered lattice data. 
These results apply to both infinite and finite spatial volume. The most general expression is 
that for a PQ theory with three sea quarks in which all valence and sea masses are different. 

This is the first calculation of a hadronic weak matrix element that is not a conserved 
current in staggered chiral perturbation theory. It illustrates the generic drawback of stag- 
gered fermions - large numbers of operators, both on the lattice and in the chiral effective 
theory, due to the additional taste degree-of-freedom. For example, while B K at 1-loop in 
continuum chiral perturbation theory is proportional to a single operator coefficient, the 
same quantity in staggered ;\PT comes from 13 operators, many of which have more than 
one distinguishable coefficient. To address this concern we show explicitly how to determine 
6 of these coefficients (at a single lattice spacing) from other simple matrix elements of the 
same lattice operator, thereby reducing the number of undetermined fit parameters. Fur- 
thermore, we suggest two ways to reduce the large number of coefficients - use of improved 
links, e.g. HYP or Asqtad, and of nonperturbative matching coefficients. As usual, partial 
quenching helps to provide more handles in the fits, and will likely be essential here. 

Use of our expression for at NLO in S^PT, in combination with sufficient lattice 
data, should allow a precise determination of B^ with staggered quarks. 

25 Recall that, generically, use of the matching procedure outlined in Eq. 1)79(1 to a given order in perturbation 
theory requires knowledge of the anomalous dimension at one higher-order. For example, 1-loop matching 
requires the 2-loop anomalous dimension. Thus, while 2-loop matching would greatly reduce the number 
of coefficients in the NLO S%PT fit, the 3-loop anomalous dimension must also be calculated. 
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APPENDIX A: B K IN CONTINUUM PQyPT WITH 2+1 SEA QUARKS 

The continuum partially quenched result for for N sea degenerate sea quarks has been 



given by Golterman and Leung 22], and checked for N sea = 2 by Becirevic and Villadoro 
The result for "2 + 1" flavors of sea quarks does not, however, exist in the literature. It can 
be obtained by taking the appropriate limit of our staggered PQ expression, Eq. (Jo"Tj). and 
here we present the result. 

To go from Eq. (JoTj) to the continuum theory, we first set all coefficients associated with 
lattice errors, i.e. those multiplied by factors of a 2 , a 2 a , a/Atr, and a 2 , equal to zero. This 
eliminates all operator coefficients but . It also removes the splittings among the sixteen 
tastes of PGBs; we are therefore free to remove the taste subscripts. Next we take the ratio 
of Bk to Bq so that all factors of four due to traces in taste space drop out. Finally we set 
m u = md ^ m s . The next-to-leading order result for Bx is then 

PQ, 2+1 



B / A8TT 2 f 2 m 2 x 



honn + hisc + bm xy + c(m 2 x - my) 2 + dm 2 xy (2m 2 D + m\ 



(Al) 

where / ~ 130 MeV. The coefficients b, c and d are unknown dimensionless low-energy 
constants proportional to the constants B, C and D in Eq. (|51|). Note that the d term is 
proportional to the sum of the sea quark masses, and, to this order, effectively renormalizes 
B in a fit at a given sea quark mass. Note also that c is absent for degenerate valence 
quarks, m x = m y . 

The chiral logarithms are given by I CO nn and Idisc- We have broken them into two contri- 
butions since the latter vanishes when m x = m y and is also the only place where sea quark 
masses enter. The result from the quark connected diagrams is 

honn = Qm%I(ml y ) - 3£(m 2 x )(m 2 xy + m\) - 3£(m 2 Y )(m 2 xy + m 2 Y ) , (A2) 
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while that for the diagrams involving a hairpin vertex is 



Idisc = I x + Iy + I v , (A3) 
(m 2 xy + m 2 x )(ml - m x ){m 2 s - m\) 

,2 



Ix = i(m x )- 

(m 2 xy + m 2 x )(m 2 D -m 2 x )(m 2 -m 2 x ) 2(m 2 xy + m x )(m 2 D - m x )(m 2 s - m\) 



[m* - m Xi 



l{m 



x, 



+ 



(m 2 — m 2 x ) 2 (m Y — m x )(m 2 — m 2 x ) 

(m 2 D - m 2 x ){m 2 - m\) - (m 2 xy + m 2 x )(m 2 - m 2 x ) - {m 2 xy + m 2 x )(m 2 D - m\) 



{ml - m\) 



,(A4) 



I Y = I X {X^Y), (A5) 

) 2 (m 2 xy + m 2 )(m 2 D -m 2 ] 
{m 2 x — mV) 2 (mY — ml)' 



2a [m\ - m 2 Y ) 2 (m 2 xy + m\){m 2 D - m 2 )(m 2 - m 2 ) 

l v — l{m ) —2 o\2f^2 • y A ®) 



The expression appears singular, but in fact is not - the poles all cancel when m 2 x , m Y — > m 2 
or m\ m\. Nevertheless, we have not found a simpler expression. We have also checked 
that this expression becomes equal to the two sea quark result (as given most compactly 



in Becirevic and Villadoro |23J) in the limit that m\ 2> m 2 D ^m 2 x ^mY^m 2 xy (so that m 2 
2m|/3). 



APPENDIX B: REVIEW OF STAGGERED CHIRAL PERTURBATION THE- 
ORY 

Here we recall the leading-order effective staggered chiral Lagrangian, whose construction 
is discussed in much greater detail in Ref . , and discuss some of its general consequences. 
The next-to-leading order Lagrangian, determined in Ref. jiol ]. turns out not to be necessary 
for B K at NLO. 

For n staggered quark flavors, spontaneous breakdown of the approximate SU (4n) chiral 
symmetry by the vacuum, 

SU(4n) L x SU(4n) R -> SU{An) v , (Bl) 
leads to 16n 2 — 1 pseudo-Goldstone bosons (PGBs). We collect them into an SU(4n) matrix, 

S = exp(i$//) , (B2) 
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where $ is a traceless An x An matrix with 4x4 submatrices: 



' U 7T+ K+ ••• ^ 



7T- D if ••• 



(B3) 



16 



Z7 = J2U a T a , etc. 



(B4) 



a=l 



and / is normalized such that f n ~ 132 Me V. We use Euclidean gamma matrices as the 
SU{A) generators: 

T a = {£5, itps, i£nv, £ M , 6} , (B5) 

where is the 4x4 identity matrix. The quark mass matrix is also of size An x An, but 
has trivial (singlet) taste structure: 



M 



mj •■• 
m d I ■•■ 
mj ■■■ 



\ 



V : 



(B6) 



/ 



Under chiral symmetry transformations, £ transforms as 

L g S't/(4n) iy , i? G S'C/(4n) Ji . 
The standard S^PT power-counting scheme is 26 



(B7) 
(B8) 



2 2 

p~ m « a 



(B9) 



so the lowest order Lagrangian is of 0(p 2 ,m,a 2 ): 



f 2 1 ?rr? 2 

^Tr(^S^St) - -uf 2 Tv{MY> + TW^t) + J£o + Dj + 5/) 2 + G 2 y ? (mo) 
o 4 o 



where Tr is the full 4n x 4n trace in both flavor and taste space, fi is a dimensionful constant 
of 0(Aqcd), and V is the taste symmetry breaking potential. 



2i i 



Note that standard S%PT does not include a in its power-counting because there arc no external operators. 
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The taste breaking potential leads to important properties of staggered mesons. For 
discussion of these features it is useful to separate V into single and double trace operators: 



V = U + W 



(Bll) 



-U 



+C 3 ^[Tr(e ) Se ) S) + / i .( 

V 

+^4^[Tr(^ ) E^E) + / i , 

V 

+C 6 ^Tr^S^St) 



(B12) 



-W = C 2y ^[Tr(^ ) S)Tr(^ ) S) + / i .c. 

V 

+^4E[ Tr ^S ) S)Tr(^ ) S) + /.. 

V 

+C 5y ^[Tr(e ) S)Tr(^ ) St)] 

V 

+C , 5 4E[ Tr (^ ) S)Tr(^ ) St)], 



(B13) 



where £^ is a 4n x 4n matrix with an ordinary 4x4 taste matrix, repeated along the 



diagonal: 



/ 



f (n) 
ST 



£ T 
£ T 

i T 



V 



(B14) 



/ 



The single trace operators in U split the tree-level PGB masses into degenerate groups: 



{ml) = (i(mt + rrij) + a 2 A F 



(B15) 



where Ap is different for each of the five S'O (4)-taste irreps P, A, T, V, and I. 27 The 
double trace operators in W generate hairpin diagrams (quark disconnected contractions) 
for flavor-neutral, taste vector and axial vector PGBs. Both U and W produce interaction 



27 



We assume in this expression that the quark flavors i and j are different. 
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vertices among even numbers of PGBs. In addition, the term, present because only the 
overall singlet can be integrated out of the theory, generates flavor-neutral, taste-singlet 
hairpins. 

It is easy to generalize the standard staggered chiral Lagrangian to quenched and partially 
quenched theories. For a PQ theory with N = n va i + n sea quarks and M = n va i ghosts, the 
chiral symmetry group is SU(AN\AM) L x SU(AN\AM) R and breaks to SU(AN\AM) V . The 
net result is that there are more PGBs, and all of the traces become supertraces (Tr — > Str) 
in the chiral Lagrangian. 

APPENDIX C: DETERMINATION OF CHIRAL B K OPERATORS 

Here we determine all of the mesonic operators which contribute to A4\ a t, and conse- 
quently to Bk-i at LO and NLO in our power-counting. There are quite a few such op- 
erators because the continuum SU(3)fi avor is enlarged to the graded, staggered symmetry 
SU(4N va i + 4iV sea |4iV m i) on the lattice. Fortunately, the fact that Ok is invariant under 
the lattice U(1)a symmetry protects the number of operators that contribute to Bk from 
becoming unwieldy. 

We first enumerate the operators which contribute to Bk at tree-level. These include the 
"standard" Bk chiral operator, and are relatively easy to determine, even in the staggered 
theory. Next we enumerate the operators which contribute to Bk at 1-loop. These can 
come either from operator-mixing on the lattice or from insertions (at the quark-level) of 
the staggered action with the Bk operator. We then determine the analytic contributions 
to Bk at NLO, relying heavily on constraints from U(1)a and CPS symmetries. Finally we 
discuss the parametric dependence of the various operator coefficients on the lattice spacing, 
a, and the strong coupling-constant, a. 

1. Operators that Contribute at Lowest Order 

We must map the staggered quark-level operator, Q^ss^red^ on ^ Q operators in the chiral 
effective theory. This is done most simply by first splitting the operator into spin vector and 
spin axial-vector parts, mapping these operators onto mesonic operators, and then taking 
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appropriate linear combinations to get chiral operators which correspond to the staggered 
Bjc operator. 

We first separate Q^ a " ere j n t two pieces 28 : 

O v = ^l(j^^)dl][^2(j^^)d2] (CI) 
O a = [^1(7^75 ® Cs)dl] [^2(7^75 ® &)d2], (C2) 

The Bk matrix element is proportional to their sum: 

Mf" ered oc< Kl\(O v + O a )\K° > . (C3) 

As in Ref. 2J|, we label these four-fermion operators by the spin and taste in each bilinear, 
such that Ov oc [V x P] and Oa oc [Ax?]. Note the convention that A = 7^5 cg> 75^, where 
the pair of spin (or taste) matrices are in separate bilinears, such that 

Of" ered oc [V x P] - [A x P] (C4) 

in this condensed notation. 

We now use spurion analysis as in Refs. Q and Q to determine the resulting chiral 
operators. The chiral structure of four-fermion operators with spins V or A is 

O f = ± (Q^(7m ® ± QZ(7m ® ^l)Ql) 2 , (C5) 

where the upper and lower signs correspond to V and A, respectively, and F^r are Hermitian 
taste matrices. When we promote the taste matrices to spurion fields, F^ r must transform 
in the following way so that Of is invariant under chiral transformations: 

F L -> LF L L\ F R -> RF R R) . (C6) 

We will set i 7 ^ and to the appropriate matrices at the end. 

We first note that chiral operators which result from two left-handed (or two right-handed) 
taste spurions have a relative minus sign when they come from a spin V versus a spin A four- 
fermion operator, whereas those that result from one left-handed and one right-handed taste 
spurion do not. Thus only operators built out of two spurions with the same handedness 
remain after taking the linear spin combination (V — A) which corresponds to Q s ^99 ered _ 



2S 



From now on we suppress color indices because both color contractions lead to the same chiral operators. 
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Operator 



±Str(Zd^ F L Zd^F L ) ± Str(EtO M SF jR St^SF jR ) 
±Str(Sa M StF L )Str(Sa M StF L ) ± StT(E%EF R )StT(E%EF R ) 

TABLE I: The two linearly-independent 0(p 2 ) operators corresponding to [(V — A) X F]. The ± 
signs correspond to spins V and A, respectively. Derivatives act only on the object immediately 
to the right. 

There are no such operators without derivatives, and two such linearly independent ones 
with two derivatives. These are shown in Table HI 

The expression in Eq. (jC5|) can be made to have the same flavor and taste structure as 
Qstaggered ^ se ^j n g ^ ne ^ wo S p Ur i ns, which we now call Fx and F 2 , equal to sparse matrices 
with a single £5 in either the sldl or sldl location. With this replacement the first operator 
in Table IJ becomes 

°x = (str(S9 At S t F 1 S9 M S t F 2 ) + Str(S t ^SF 1 S t ^ t SF 2 )) (C7) 
n 

and the second becomes 

°x = H (str(S9 At StF 1 )Str(S9 At StF 2 ) + Str(St9 ft SF 1 )Str(St^SF 2 )) . (C8) 

bmce Of has the same left-left current structure as the standard continuum full QCD 
operator ^J, it should contribute to Bk at tree-level, and it does. The matrix Fx annihilates 
the Kx (an sldl meson) while _F 2 produces the K\ (a d2s2 meson). 

Because Bk in continuum xPT only depends on a single parameter the coefficients 
of and cannot be independent. We can determine their precise relationship using 
the fact that both operators in the chiral effective theory arise from the same quark-level 

• /rtstaqqered 

operator, U K yy : 

^taggered SxPT^ ^ qK + ^ Q ^ (cg) 

but they have nonvanishing tree-level matrix elements between different types of kaons. In 
particular, while only has a tree-level matrix element between K® and (by construc- 
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tion), only has one between "mixed" kaons, i.e. sld2 and s2dl: 

<^w+wo = -^- (cm 

By comparing the same two matrix elements at the quark level we can extract the ratio 
C£ : . As written in Eq. Q s ^ a " ered can contract with K® and K®, but not the mixed 
kaons. However, it can be Fierz-transformed into an operator with the appropriate mixed 
flavor structure, at the cost of introducing new tastes. Using the Fierz-transformation rules 
in Appendix A of Ref . 2(j] , we find that 

[(V -A)xP] ^> - A) x (S + P + T - V - A)]. (C12) 

Thus the ratio of the two matrix elements is 

(Kf p \Of" ered \K%) 



/jy-Q"""'\/nstaggered\ j^Qmix\ 
\ A 1P \ W K \ JX 2P / 



(C13) 



and = . The Fierz transformation also reveals why the seemingly "extra" operator, 



, is necessary in the chiral effective theory - it accounts for the fact that Q^sa^ed a 
nonvanishing tree-level matrix element between flavor-mixed kaons of all tastes. 



2. Operators that Contribute at NLO 

Bk receives NLO contributions from a host of additional operators because of lattice 
discretization effects. In order to determine all such operators we must consider the three 
different ways in which they can arise at the quark-level: through perturbative matching, 
mixing with higher- dimension operators, and insertions of operators from the staggered 
action. 



a. 1-Loop Matching Coefficients 

Lattice operators can mix with each other as long as they are in the same representation 
of the symmetry group that maps a hypercube onto itself. Consequently, lattice operators 
which correspond to continuum operators with a particular spin and taste mix with other 
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operators that correspond to a different continuum spin and taste. Specifically, although 
Qstaggered ^ ^ e continuum is strictly a [{V — A) x P] four-fermion operator, once it is 
put on the lattice it mixes with [(V + A) x P] and other spin-taste structures. Like all 
taste-breaking effects, this operator mixing is due to hard gluon exchange, so the amount 
of mixing can be calculated perturbatively in a. Such 1-loop "matching coefficients" for 
staggered four-fermion operators have been calculated for both naive fermions and using 
hypercubic fat links Q]. Some of this operator mixing is accounted for and removed in 
present lattice calculations of Bk ■ It is now standard to use 1-loop matching coefficients 
to remove the (9(a/47r) mixing between the desired [(V — A) x P] lattice operator and the 
unwelcome [(V+A) x P] lattice operator. This still leaves mixing at 0(a 2 ) with [(V+A) x P] 
and at 0(a/4ir) with wrong taste operators. We discuss them in turn. 

Lattice operators with spin structure (V + A) contribute to Bk , albeit suppressed by 
0(a 2 ). Recall from the previous section that such operators arise from the combination of 
one left-handed spurion and one right-handed spurion, where Fl and Fr transform as in 
Eq. l)C6j) . Here we need only consider operators without derivatives. This is because any 
two-derivative operators are of 0(a 2 p 2 ) and can only contribute to Bk at NLO through 
analytic terms, which we determine separately in Appendix IC 31 There is only one operator 
without derivatives, 

O 1 / = Str(F 1 SF 2 St), (C14) 

and it cannot contribute to Bk until 1-loop because of its flavor structure. Note that the 
1-loop contribution of 0^, p is of NLO in our power-counting because its coefficient is of 
0(a 2 ) and the loop momentum further suppresses it by 0(p 2 ). 

Lattice calculations to date have not taken into account the 0(a/4ir) mixing between 
\{V — A) x P] and lattice operators of different tastes. Bk therefore receives contributions 
from "wrong taste" lattice operators, but suppressed by 0(a/47r). Using the matching 
coefficients from Ref. Q| we find that the desired lattice operator mixes with the following 
eight four-fermion operators at 1-loop: 

[VxS], [AxS], [VxT], [AxT], [SxV], [SxA], [PxV], [PxA], [TxV], [TxA], 

(C15) 

as well as taste off-diagonal operators, i.e. those in which the two bilinears have different 
tastes, and operators which violate rotational and S'O (4)-taste symmetry. We map all of 
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Four-fermion Op. Chiral Op. 

[V x S] and [AxS] — ► Str^/E^E*) 
[V x T] and [A x T] — > Str(^E^E+) 
[5 x V] and [P x V] — ► Str(^E)Str(^Et) 

Str(^E)Str(^E) + Str(^Et)Str(^Et) 

Str(^E£ M E) + Str(^Et^Et) 
[S x A] and [PxA] — > Str(^ 5 E)Str(^Et) 

Str(^ 5 E)Str(e 5M E) + StrfotfE^Str^Et) 

Str(e M5 Ee 5M E) + Str(e M5 StC 5M Et) 

TABLE II: Mesonic operators corresponding to the four-fermion operators listed in Eq. IC151 Re- 
peated indices are summed with the constraint that jjL 7^ v in the taste tensor matrices. 

these quark-level operators onto chiral operators using a straightforward generalization of 



the spurion analysis in Ref. 24| to flavor off-diagonal operators. TablelTDcontains the results. 
We do not show the chiral operators generated by [T x V, A] because they are simply a subset 
of those from [S, P x V, A] . Nor do we include taste off-diagonal operators, as they do not 
contribute to any PGB quantities of interest. Furthermore, we neglect operators which break 
rotational and/or 50(4)-taste symmetry, as doing so would require at least four repeated 
Lorentz indices in the chiral operator. This would have to come from the presence of four 
derivatives, four additional taste matrices, or two derivatives and two taste matrices, and 
therefore be of a higher order than we consider here. 

Because these operators arise from mixing with Q s ^ a 99ered ^ faey have the same flavor 
structure as the original continuum operator. Thus ^ — > Fit and £ v(i — > F 2 t, where Fit 
(F 2 t) is a sparse matrix with a single ^ (£1^) in the sldl (s2d2) location. We define Fu, 
Fu, Fi V , F 2V , Fia, and F 2 a analogously, such that the following eight operators contribute 
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to B K : 29 



0" 


= J]Str(F i; EF 2 ,Et) 


(C16) 


Of 


= ^Str(F 1T EF 2T E t ) 


(C17) 


Of 


= ^ (str(F iy £F 2y £) + Str(F iy S t F 2K S t )) 


(C18) 


of 


= ^ (str(F iy £)Str(F 2y £) + Str(F iy E t )Str(F 2V E t )) 

/V, 


(C19) 


of 


= ^Str(F iy S)Str(F 2y E t ) 


(C20) 


of 


= ^ (str(F 1A £F 2A £) + Str(F 1A StF 2A St)) 


(C21) 


of 


= (str(F 1A £)Str(F 2A £) + Str(F 1A St)Str(F 2A St)) 


(C22) 


of 


= ^Str(F 1A £)Str(F 2A £t). 


(C23) 



Moreover, since Eq. (j(J15|) turns out to be a complete list of U(1)a and S'0(4)-rotation 
invariant four-fermion operators, these are also the only wrong-taste operators that can 
contribute to at 0(a 2 ). Thus all of the operators in Eq. (jC23|) have two independent 
coefficients, one of 0(a/4ir) and the other of 0(a 2 ). Note that, as one would expect, such 
"wrong taste" operators can only contribute to the Bk matrix element through loops. Such 
contributions are of NLO in our power-counting because of the factor of a/47r or a 2 in the 
coefficients and the further suppression of 0(p 2 ) by the loop- momentum. We therefore do 
not need to consider wrong-taste operators with derivatives. 



b. Dimension 8 Quark-level Operators 

The four-fermion operator Q s ^ L 99 ered j s dimension 6. As discussed in the previous sub- 
section, it mixes with other dimension 6 operators on the lattice, but such mixing is due to 
gluon exchange and suppressed by a/4n. Therefore the naively discretized version of the 
continuum operator, Q s £ a 99 ered < lat ; generates the dominant contribution to the lattice K^—K® 



29 We note as an interesting aside that O^ 1 does not appear in the staggered chiral Lagrangian because, if 
£r is flavor diagonal, it is just a trivial constant. 
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matrix element, Aiiat- However, many higher-dimension quark-level operators respect the 
same lattice symmetries as Q s ^ a 99 ered anc [ thus can also contribute to Bk on the lattice. 
Dimension 7 and 8 quark-level operators are explicitly suppressed in the action relative to 
dimension 6 operators by factors of a and a 2 , respectively, so they can map onto chiral oper- 
ators which generate 1-loop contributions to Bk at NLO. In fact, because a 2 ~ a 2 ~ a/47r 
in our power-counting, they are just as or more important than the dimension 6 operators 
which contribute through mixing. Fortunately we need not go higher than dimension 8 at 
the quark-level because 1-loop contributions to Bk from dimension 9 operators are at least 
0(a 3 p 2 ) and thus of higher than NLO. 

At the quark level, dimension 7 operators contain four fermions plus either a derivative or 
a factor of the quark mass matrix. However, all dimension 7 quark level operators turn out to 

any PGB processes of interest. 

Thus we move on to dimension 8 quark-level operators, which contain four fermions plus two 
derivatives, two factors of the quark mass matrix, or one of each. Dimension 8 operators with 
one derivative and one mass matrix must have a single nontrivial spin or taste matrix (either 
V or A) in order to contract up the free derivative index. Such operators are either spin 
off-diagonal or taste off-diagonal; both are irrelevant. Dimension 8 operators with two mass 
matrices map onto chiral operators with at least two mass spurions. Because they also carry 
a factor of a 2 from the action, such operators are of 0(a 2 m 2 ) or higher in our power-counting 
and can be neglected. We therefore need only consider dimension 8 four-fermion operators 
with two derivatives that possess the same lattice symmetries as Q s ^ a 99 ered _ Enumerating 
all such operators would be a tedious exercise in staggered group theory. Fortunately, this 
exercise is unnecessary if we start with all chiral operators which could contribute to Bk at 
NLO and work backwards. 

All chiral operators which arise from dimension 8 quark-level operators automatically 
come with a power of a 2 . Thus, in order to generate 1-loop contributions to Bk at NLO, 
they cannot contain any derivatives or be suppressed by additional powers of a or a. All of 
the operators listed in Table ITT1 satisfy this condition. In fact, it turns out that there are no 
additional chiral operators which are invariant under U(1)a and derivative- free. Therefore all 
of the chiral operators in Eqs. (|C14|) and (|C16|) - ()C23|) should contribute to Bk both through 
1-loop mixing, suppressed by C(a/47r), and through dimension 8 operators, suppressed by 
0(a 2 ). 
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Yet there is a subtlety - not all a 2, s are created equal. This is because, at the quark level, 
taste-changing four-fermion interactions require gluon exchange, and therefore have hidden 
powers of a which do not appear explicitly in the chiral effective Lagrangian, Eq. (jBlOj) . We 
must therefore consider the 0( u a 2 ") operators in more detail to allow correct extrapolation 
of lattice data. 

In the case of matching the continuum version of Q s £ a 99 ered on ^ dimension 8 lattice 
four-fermion operators, we want to find the dimension 8 operators that match onto with 
\{V — A) x P] without gluon exchange, because they truly contribute at 0(a 2 ). Any other 
dimension 8 operators contribute at 0(aa 2 ) or higher, and thus are of at least NNLO in our 

n o 

power counting. As discussed in both Refs. |29|] and [30(, lattice bilinears of different spins 
and tastes "mix" 30 with each other at tree-level in perturbation theory, but to all orders in 
a, simply because of the fact that the spinor components of staggered fermions are spread 
out over a hypercube. 31 Fortunately, at tree-level, one can consider the two bilinears in each 
four-fermion operator separately. Thus, in order to determine which four-fermion operators 
mix at 0(a 2 ), one needs only find which bilinears mix at 0(1), 0(a), and 0(a 2 ). The 
product of two 0(a) bilinears, or of one 0(1) and one 0(a 2 ) bilinear, will be a four-fermion 
operator that mixes at O (a 2 ). 

n 

We do so using Eqs. (26) and (27) in Ref. [30j|. We find that the dimension 3 bilinear 
[7^(1 + 75) ® £5] mixes with taste V dimension 4 bilinears at 0(a). It also mixes with 
dimension 3 and 5 taste P bilinears at 0(1) and 0(a 2 ), respectively. Because spin off- 
diagonal four-fermion operators cannot contribute to PGB processes, we are restricted to 
multiplication of each of these bilinears by itself. For example, we can multiply the bilinear 
[1 x £ M ] times itself to generate the four-fermion operator 



The above operator is simply a generalization of the dimension 6 four-fermion operator 
[S x V] to dimension 8. The derivatives give it the appropriate dimension, and allow it to 
have both SO (4) rotational symmetry breaking (/i = v) and conserving (fi 7^ v) pieces. 

30 Recall that we are using the terminology "mix" as shorthand for this process of matching lattice operators 
onto those in the continuum. 

31 Note that these mixings are independent of both the choice of gauge links and the choice of how to 
construct gauge- invariant operators. 




(C24) 
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>cx-x 



FIG. 3: Schematic of B^ vertex modification due to the insertion of an operator from the staggered 
action. In the first diagram the square represents the original Bk vertex and the circle represents 
the inserted strong four fermion vertex. In the second diagram the square represents the resulting 
local vertex in the effective theory. 

It is clear now that the lattice version of Q^ a 99 ere really only mixes at 0(a 2 ) with 
dimension 8 quark- level operators that have spin-taste structures [S x V] , [P x V] , [TxF], 
\{V — A) x P], and [{V + A) x P] along with two derivatives. These lead to the four chiral 
operators Oi v , 2V , Ol v , and 0}, p . 32 We therefore conclude that their coefficients are of 

X X X X 

0(a 2 ), while the rest of the operators enumerated this section can be neglected. This is 
summarized in the second column of Table IIV1 



c. Insertions of the Staggered Action 

The lattice operator mixing discussed in the past two subsections is not particularly 
intuitive because it is a stems from perturbation theory and complicated lattice symme- 
tries. In contrast, the fact that operators from the staggered action modify the quark-level 
P>k operator, Q s ^ a 99 ered ^ can De snown through a simple diagrammatic argument. Consider 
the four-fermion Q s ^ a 99 ered anc i another four-fermion operator from the staggered action. 
Two of the quarks in Q s ^99 ered can contract with two of the quarks in the other operator to 
form a 1-loop diagram with four external quarks, as shown in Fig. 01 This produces a local 
operator in the chiral effective theory because the loop propagates over 0(1/Aq C d), rather 
than over 0(l/m n ). It corrects Q s ^ a 99 ered a ^ 0(0,^) because all of the dimension 6 four- 
fermion operators in the staggered action are taste-breaking. Such 0{a 2 a ) corrections turn 
out to be numerically enhanced, however, so they are still of NLO in our power-counting. 
We must therefore determine all chiral operators produced by combining Q s ^ a 99 ered wr th each 



32 They also lead to higher-order chiral operators which break 5*0(4) taste symmetry, corresponding to the 
piece of Eq. (|C24f) with u = v. 
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of the dimension 6 four-fermion operators in the staggered action. 

Ref. Q| developed a method for combining four-fermion operators at the chiral level. 
Essentially, one forms chirally invariant operators which contain two taste spurions from the 
first four-fermion operator and two from the second. Q s ^ a 99 ered spin-taste structure [{V — 
A) x P], so the two taste spurions from Q s ^9s ered m ust transform as in Eq. (jC6|) . Moreover, 
we learned in Sec. IC II that only operators with two left-handed (or two right-handed) taste 
spurions survive after taking the linear spin combination (V — A). In contrast, the four- 
fermion operator may have any spin - V, A, S, P, or T - combined with an appropriate taste 
to maintain overall £/(l)A-invariance. Furthermore, the two taste spurions from the inserted 
four-fermion operator must be flavor diagonal because the action does not change quark 
flavor. 

Operators with spins V or A can have tastes S, P, or T, and their taste spurions transform 
as Eq. (|C6|) . Generically, there are three possible ways of combining these spurions with 
those from U K : 



(F L F L )(F> L F' R ) + p.c. 
{F L F L ){F' L F' L ) + p.c. 

(F L F L )(F> R F R ) + p.c., (C25) 

where "p.c." indicates parity-conjugate. Here the unprimed spurions are from Q s ^ a 99 ered 
and the primed are from the four-fermion operator in the staggered action. Because the 
two Fl's will ultimately have different off-diagonal flavor structures, they must be separated 
by at least one E-field on either side within a supertrace. We use this fact to eliminate a 
few operators. Consequently, these spurion combinations only lead to two nontrivial generic 
chiral structures, shown at the top Table 11111 To produce operators which contribute to 
B K we let F L — > Fi$zF 2 and F' R — > £j,£s, or It is easy to see that the chiral operators 
are trivial when the taste spurions from the quark-level operator in the action are just the 
identity. Thus we are left with four new operators which contribute to B>k with coefficients 
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Generic Structure 


Specific Operator 


StT(F L XF^F L XF' R tf) + p.c. 


-> Str(FiF 2 ) +p.c. 




Str(e 5 SF 1 Ste 5 SF 2 St)+p. c . 




Str(^EF 1 Et^S J F 2 Et)+ p.c. 


Str(F L EF^Et)Str(F L EF^Et) + p.c. - 


Str(Fi)Str(F 2 ) + p.c. 




Str(e 5 SF 1 Et)StrfeEF 2 Et) +p. c . 




Str(^SF 1 St)Str(e 5 S^St) + p. c . 


§tT{F L F L tfF L Y,F R )+p.c. 


- Str(F 1 ^StF 2 S^)+p.c. 




Str(F 1 e At5 St J F 2 E^ M )+p.c. 


StT(F L F L ^)StT(F L EFR) +p.c. 


- Str(F 1 C M Et)Str(F 2 E^)+p.c. 




Str(F 1 ^ 5 Et)Str(F 2 E^)+p.c. 


Str(F L F L StF i ?iSt)+p.c. 


- Str(F 1 ^EtF 2 ^Et)+p. c . 




Str(F 1 ^ 5 EtF 2 ^Et)+p. c . 


Str(F L F L St)Str(F L F^St) + p.c. 


- Str(F 1 e At Et)Str(F 2 e M Et) +p. c . 




Str(F 1 ^ 5 Et)Str(F 2 ^Et) +p. c . 



TABLE III: Mesonic operators corresponding to insertions of four-fermion operators from the 
staggered action. F\ and F% are flavor off-diagonal and come from the Bk quark-level operator. 
The remaining taste matrices are flavor-diagonal. Repeated indices are summed with the constraint 
that fi ^ v in the taste tensor matrices. The notation "p.c." indicates the parity-conjugate of the 
previous operator. The difference between upper and lower panels is described in the text. 

ofO(a 2 J: 

Of = Str(e 5 EF 1 S t e5SF 2 S t )+p.c. (C26) 
Of = Str(e 5 EF 1 E t )Str(e5EF 2 E t )+p.c. (C27) 

°x = Yl Str(^SF 1 St)Str(^SF 2 St) + p.c. . (C29) 
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Operators with spins S or P can have tastes V or A 33 Their chiral structure is as follows: 

0' F = (qI(1<8> F l )Qr ± Qr~(1 ®Fr)Ql) 2 , (C30) 

where and the upper and lower signs correspond to spin S and P, respectively, so their taste 
spurions must transform as 

F L -> LF L R\ Fr - RFrL ] . (C31) 
As before, there are three possible ways of combining these spurions with those from 

sr\staggered 

u K ■ : 

(F l F l )(F l Fr) + p.c. 
(F L F L )(F L F r L )+p.c. 

(F L F L )(F R F R )+p.c. (C32) 

These spurion combinations generate the four nontrivial chiral structures shown in the lower 
panel of Table IIHI In this case, to produce operators which contribute to E>k we let Fl — > 
Fi$zF 2 and F L , F R — > £ M or Because two of the taste matrices are diagonal, the following 
(and similar) relations allow further simplification of the resulting operators: 





= ~F 1 




= F 1A 




= F 1V . 



(C33) 

Such simplification reveals that these operators are none other than other than 0\ p and O 1 ^ 
- 3A . However, here they have arisen through insertions of the staggered action rather than 
lattice operator mixing. 

In principle, we must also consider insertions of the dimension 4 quark mass term, which 
can modify Q s ^ a 99 ered j n a s j m ilar manner to the four-fermion operators already considered. 
The resulting chiral operators would contain the two taste spurions from Q s ^ a 99^ed a 
mass spurion, and thus be of 0(m) in our power-counting. However, it turns out that there 
are no such operators for quite a simple reason. Recall that, since F\ and F 2 are flavor 



33 We choose to neglect spin T in this discussion because it ultimately leads to a subset of operators already 
generated by spins S and P. 
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off-diagonal, they must either be separated by a S-field on both sides or be in different 
supertraces, each with a E. A single mass spurion is insufficient for meeting this criterion. 
Therefore operators with mass spurions only contribute to Bk at NLO through analytic 
terms, which we discuss next. 



3. Analytic NLO Contributions 

We finally consider analytic NLO contributions to the Bk matrix element, M.k- Many 
operators generate contributions with the same quark mass dependence; it is unnecessary to 
separate them in fits to lattice data. Thus rather than enumerating all possible operators, we 
can determine all linearly-independent functions of the quark masses, allowing symmetries 
to do much of the work. 

Given our power-counting scheme, generically six types of operators can contribute to 

0(p 4 ), 0("a 2 V), 0("a 4 "), £>(m 2 ), 0(p 2 m), £>("a 2 "m), (C34) 

where quotation marks indicate that a 2 can be interchanged with either a 2 , a/4ir, or a 2 a . 
We can immediately rule out contributions from 0(p A ) operators because there are four 
derivatives but only two fields upon which they can act. We can also rule out NLO ana- 
lytic contributions from (9("a 4 ") operators. Without derivatives or mass matrices the only 
possible matrix element from such operators is cx "a 4 " . However, this contribution does not 
vanish in the chiral limit, and therefore violates the U(1)a symmetry possessed by Q s ^ a 99 ered ^ 
so it cannot occur. Effectively, U(1)a forces all tree-level matrix elements between taste-^5 
kaons to vanish in the chiral limit. 34 

Determination of the NLO analytic contributions to Aix from (9("a 2 "p 2 ) operators is 
also straightforward. The derivatives must act on the two external kaons, bringing down 
a factor of p 2 K which becomes m 2 K when the kaons are on-shell. Therefore the first NLO 
analytic term in the expression for Aix is 

Am 2 K , (C35) 



34 Of course, this is nothing more than a restatement of the fact that taste-^5 mesons are lattice Goldstone 
bosons, but we have nevertheless checked this result explicitly. 
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where A is an undetermined coefficient that is, in principle, suppressed by either a 2 , a 2 , 
a/ '47T, or a 2 a . For fitting it is important to know whether or not A does, in fact, depend upon 
all four expansion parameters. Operators of 0{a 2 p 2 ) come from mixing with dimension 8 
operators. Recall from Appendix IC 2 bl that [(V — A) x P] mixes with taste P and taste V 
operators at 0(a 2 ). The operator itself comes from a taste P quark-level operator, so 
its coefficient, C£ , receives a correction of 0(a 2 ). One can also easily map the quark-level 
operator [(V + A) x P] onto a two-derivative chiral operator that contributes to KA K at 
tree-level: 

Str(S^S t F 1 )Str(S t 9 M SF 2 ). (C36) 

Although we lump all NLO analytic contributions to M.k of the form a 2 m 2 K together, we note 
that it would be important to separate the correction to from the rest if one were to do 
a fit including some NNLO terms. Operators of 0(a 2 p 2 ) come from perturbative matching 
with all tastes of four-fermion operators on the lattice, which we discussed in Appendix IC 2 al 
As we just showed, mixing with taste P four-fermion operators leads both to corrections to 
C£ and to a new operator. This time they produce contributions to M. K of the form a 2 m 2 K . 
In contrast, operators of 0(p 2 a/4ir) only come from perturbative matching with wrong- 
taste dimension 6 four-fermion operators. Because wrong-taste operators cannot contribute 
to M.k at tree-level, however, there is no analytic contribution of the form m 2 K a/4TT. Finally, 
operators of 0(p 2 a 2 ) come from insertions with operators in the action. The insertion of 
QsUiggered w ^ ^ e taste-singlet dimension 6 four-fermion operator maps onto , but now at 
higher-order. Thus there is an NLO analytic correction to Cr; , and therefore a contribution 
to M.K-, of the form a 2 a m 2 K . We summarize the parametric dependence of the coefficient A 
in Table El 

Determining the analytic NLO contributions to Aix from the last three types of operators 
requires a bit more work. To do this, we utilize a sy mmetry of Q s ^ a 99 ered w Yn.oh has not 
been explicitly needed until now: CPS symmetry |27[. Recall that, because Q s ^ a 99 ered j g a 
weak operator, it is not parity-invariant, as is the strong chiral Lagrangian. However, it 
is invariant under CPS symmetry, and this can be used to restrict the chiral operators to 
which is corresponds. Standard CPS symmetry is CP plus d ^ s and ^ m s . In our case, 
because there are two sets of valence quarks, we must impose a modified CPS symmetry 
with both dl <-> si and d2 <-> s2 so as to make the quark-level operator invariant. At the 
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chiral level, this corresponds to the following transformations: 

C : E <-> E r 
P : E w E f 
"S" : dl «-> si, d2 w S 2, 

where the last line denotes an interchange of indices in both E and M. One can easily show 
that all of the operators enumerated in the previous subsections are already invariant under 
this series of transformations. However, CPS symmetry will prove important for restricting 
the analytic contributions to M.r- 

First consider analytic contributions to Mr from 0(m 2 ) operators, which necessarily 
contain the two spurions from Q s ^99 ered pj us two m ass spurions. Because of the flavor 
structure of Ft and F 2 , such operators will only produce nonvanishing tree-level matrix 
elements if one mass spurion is within each supertrace, e.g.: 

Str(F 1 EM t )Str(F 2 EM t ) + Str(F 1 ME t )Str(F 2 ME t ) (C37) 
Str(F 1 EM t )Str(F 2 ME t ) + Str(FiME + )Str(F 2 EM t ) . (C38) 

Note that the second term in each operator is the CPS conjugate of the first term. When 
contracted with an external F 2 and a the first operator generates an analytic term 
proportional to m d + m 2 s , while the second generates one proportional to mdTn s . We choose 
to combine these terms such that the new analytic contributions to M.r at NLO are of the 
form 

Bm%, C{m d - m s ) 2 . (C39) 

Note that the second term vanishes quadratically as (m^ — m s ) — > 0. 

Next consider 0(p 2 m) operators which consist of the two spurions from Q s ^ a " er j two 
Lie derivatives, and a mass spurion. To produce a nonvanishing matrix element, each Lie 
derivative must be in a separate supertrace with either Fi or F 2 , e.g.: 

Str(F 1 E^E t )Str(F 2 E^E t )Str(EM t + ME f ) , (C40) 
Str(F 1 E^E t )Str(F 2 E^E t EM t ) + Str(F 1 E^E t )Str(F 2 ME t E^E t ) . (C41) 

Once again, these operators are CPS-invariant. The first operator generates a term pro- 
portional to m| r Str(M) = m| r Tr(M sea ), while the second generates one proportional to 
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m l'( m rf + m s) m K- Thus there is only one new analytic contribution to M.k at NLO from 
such operators: 

Dm 2 K Tr(M sea ). (C42) 

Lastly we consider (9("a 2 "m) operators. These contain two spurions from Q s ^ a 99 ered ^ 
possibly two flavor- diagonal taste spurions, and a mass spurion. Thus their tree-level ma- 
trix elements are proportional to one power of the quark mass, for which there are three 
possibilities: {m d + m s ), {m d — m s ), and Ti(M sea ). However, "a 2 "(m rf — m s ) is not CPS- 
invariant, and "a 2 " (m^ + m s ) is included in the parameter A, so we need only consider 
analytic terms of the form u a 2 "Ti(M sea ). The factor Tr(M sea ) always arises from the op- 
erator Str(EM' + MET) wri en E = 1, so the only operators that can produce tree-level 
contributions of the form "a 2 "Tr(M sea ) are simply the previously enumerated 0(a 2 ), 0(a 2 ), 
0{a/4ir), or 0(a 2 a ) operators multiplied by Str(EAf f + MY)). However, C("a 2 ") operators 
cannot generate tree-level contributions to M.r because of the U(1)a symmetry, so neither 
do these particular (9("a 2 "m) ones. 



4. Parametric Dependence of Operator Coefficients 

Many of the chiral operators enumerated in the previous sections come from multiple 
quark-level operators. Thus they can have multiple coefficients, each multiplied by a different 
expansion parameter. While this does not affect the form of at NLO in S^PT, it 
certainly affects extrapolation of lattice data using multiple lattice spacings. Thus, for 
ease of use, we summarize the dependence of the NLO operator coefficients on the various 
expansion parameters in Table II VI All operators which come from 1-loop operator mixing 
have coefficients of 0(a 2 ), 0(a/Aiv), or both, as described in Section IC 2 al Operators which 
arise through mixing with dimension 8 four-fermion operators come in at O (a 2 ). Finally, 
operators that arise from insertions of four-fermion operators are of 0(a 2 a 2 ). Thus the 
various Bx operators have anywhere between 1 and 4 independent undetermined coefficients. 
For example, O}^ has four coefficients, multiplied by a/47r, a 2 , a 2 , and a 2 , respectively, 
whereas C 2P has a single coefficient multiplied by a 2 . 

Finally, for completeness, we include the dependence of the various analytic terms on a 
and a in Table IVl 
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Bx Operator Source - Coefficient Dependence 

1-loop Mixing? Dimension 8 Op.? Action Insertion? 



u x 


Yes — ct 2 




Yes — oP 


Vpc — rp> 
J. CO IAq, 


QlT 

^X 


Yps — rv Ait 

i. CO LX 1 /I , 


a 2 


No 


No 


^x 


Yes — c\ / 4-tt 

J. V/O '.1 I ^Cl I , 


a 2 


Yes — a, 2 


Yes — n 2 


Q2V 
^X 


YP<? — rv /4-7T 

J. CO LX 1 /I , 


a 2 


Yes — op 


Yps — 

J. CO Q 


Q3V 
^X 


Yes — /4-7T 


a 2 


Yes — cl 2 


Yps — cP 1 


QlA 

^X 


Yes — (y /4-TT 


a 2 


No 


Yes — 

ico U, Q 


U X 


Yp<2 — /'V / ZlTT 
X CO Lx / ^til ■ 


a 2 


No 


Vpc — rp 1 
-L CO CXq, 


°x 


Yes — olI Att , 




1NO 


Yes — a^. 




Yes - a/47r , 


a 2 


No 


No 




No 




No 


Yes - a? a 


^x 


No 




No 


Yes - a„ 


Q 2T 
^X 


No 




No 


Yes - a 2 a 


^X 


No 




No 


Yes - a„ 



TABLE IV: Parametric dependence of NLO Bk operator coefficients. Each chiral operator can 
come from more than one quark-level operator, and therefore enter the NLO expression for 
Bk multiplied by more than one expansion parameter. 
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Br Analytic Term Coefficient Dependence 

Am 2 K a 2 , a 2 , 

Bm% 1 

C(m d - m s ) 2 1 

Dm 2 K Tr(M sea ) 1 

TABLE V: Parametric dependence of NLO analytic contributions to Bk • An entry of "1" indicates 
that the coefficient is a constant at this order. 
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